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Executive  Summary 

Vehicular-towed  magnetometer  arrays  have  been  used  for  MEC  detection  since  the  late  1980s. 
However,  most  vehicles  are  highly  ferromagnetic  due  to  their  ferrous  frame,  skin,  and  drive  train, 
and  the  resulting  magnetic  self- signature  can  easily  overwhelm  the  signal  from  subsurface 
objects  and  render  the  data  useless.  Further,  because  the  vehicle  signature  is  induced  by  the 
Earth’s  magnetic  field,  it  is  not  constant;  it  changes  primarily  with  the  vehicle’s  orientation 
relative  to  North,  and  secondarily  with  the  vehicle’s  pitch  and  roll.  Several  successful  vehicle- 
towed  magnetometer  arrays  have  addressed  the  vehicle  signature  problem  through  the  use  of 
custom-built  non-ferrous  aluminum- framed  vehicles  which  minimize  vehicle  self-signature. 
However,  the  cost  of  these  vehicles  was  in  excess  of  $100k,  putting  them  out  of  range  of 
commercial  UXO  contractors.  The  logical  question  is:  is  this  kind  of  expensive  custom  vehicle 
absolutely  necessary  to  acquire  high-quality  towed  array  magnetometer  data,  or  can  a  contractor 
employ  a  vehicle  with  a  higher  signature  and  filter  out  its  effects?  Under  this  project  we  tested  a 
number  of  commercial-off-the-shelf  (COTS)  side-by-side  utility  vehicles  (UTVs)  and  an  all- 
terrain  vehicle  (ATV)  for  their  applicability  as  tow  vehicles  for  a  towed  magnetometer  array  by 
measuring  their  magnetic  signature  and  seeing  if  the  signature  can  be  removed  through  simple 
filtering  techniques  to  yield  data  of  a  similar  quality  to  data  obtained  using  a  custom-built 
vehicle.  SAIC  selected  what  it  felt  was  the  best  compromise  of  low  signature,  cargo  space, 
terrain- handling  capability,  and  ride  -  the  aluminum- framed  Club  Car  XRT  1550  -  purchased  it, 
and  adapted  it  to  pull  SAIC’s  full-sized  VSEMS  mag/EM61  towed  platform  (developed  under 
ESTCP  project  MM-0208).  We  found  that  this  vehicle  engenders  only  a  moderate  signature  in 
the  data,  and  that  this  signature  can  be  easily  removed  with  the  de-median  filter  in  Geosoft  Oasis 
Montaj  that  is  already  commonly  employed  for  filtering  geophysical  data. 

1  Introduction 

1.1  Background 

Total  field  magnetometers  are  proven  instruments  for  MEC  detection.  Their  high  sensitivity 
makes  them  particularly  appropriate  for  detection  of  major  caliber  air-dropped  munitions  like 
2501b  bombs  that  are  large,  ferrous,  and  can  penetrate  deep.  In  order  to  maximize  the  number  of 
survey  acres  per  day,  a  wide  swath  may  be  obtained  by  ganging  total  field  magnetometers 
together  in  an  array.  Because  of  the  weight  of  the  required  batteries  and  electronics,  the  wear  and 
tear  on  the  person  carrying  the  equipment,  and  the  desire  to  increase  survey  speed,  it  is  natural  to 
want  to  tow  such  a  magnetometer  array  behind  a  vehicle.  However,  the  magnetometer’s  very 
sensitivity  makes  this  difficult.  Most  vehicles  are  highly  ferromagnetic  due  to  their  ferrous 
frame,  skin,  and  drive  train,  and  a  towing  vehicle’s  magnetic  self- signature  can  overwhelm  the 
signal  from  subsurface  objects  and  render  the  data  useless.  Further,  because  the  vehicle  signature 
is  induced  by  the  Earth’s  magnetic  field,  it  is  not  constant;  it  changes  primarily  with  the  vehicle’s 
orientation  relative  to  North,  and  secondarily  with  the  vehicle’s  pitch  and  roll.  Historically, 
successful  vehicle-towed  magnetometer  arrays  (e.g.,  STOLS,  MTADS  and  VSEMS)  have 
addressed  the  vehicle  signature  problem  through  the  use  of  vehicles  custom-built  out  of  non- 
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ferrous  materials  to  minimize  vehicle  self-signature.  Both  MTADS  and  VSEMS  utilize  dune 
buggies  custom-built  by  Chenowth  Racing  with  an  aluminum  frame  and  a  magnesium  alloy 
engine.  Fifteen  years  ago,  the  cost  of  these  vehicles  was  in  excess  of  $100k,  putting  them  out  of 
range  of  most  commercial  UXO  contractors.  The  logical  question  is:  is  this  kind  of  vehicle 
absolutely  necessary  to  acquire  high-quality  towed  array  magnetometer  data,  or  can  a  system 
employ  a  vehicle  with  a  higher  signature  and  filter  out  its  effects? 

Since  the  development  of  MTADS  and  VSEMS,  a  variety  of  commercial-off-the-shelf  (COTS) 
all-terrain  vehicles  (ATVs)  and  so-called  “side-by-side  utility  vehicles”  (UTVs)  have  become 
available.  An  ATV  is  a  small  vehicle  that  is  straddled  like  a  motorcycle  and  has  handlebars  and 
hand  controls  like  a  motorcycle,  whereas  a  “side-by-side  UTV”  looks  like  a  cross  between  a  golf 
cart  and  a  small  Jeep,  typically  has  a  small  dump  bed  on  the  back,  has  upright  seating  for  two 
adults,  has  the  seats  oriented  side-by-side,  and  has  a  steering  wheel  and  brake  and  accelerator 
pedals  like  a  car.  By  way  of  example,  a  ubiquitous  side-by-side  UTV  is  the  John  Deere  Gator. 
The  objective  of  this  project  was  to  test  a  number  of  these  COTS  vehicles  for  their  applicability 
as  tow  vehicles  for  a  towed  magnetometer  array  by  measuring  their  magnetic  signature  and 
seeing  if  non-zero  vehicle  signature  can  be  removed  through  simple  filtering  techniques  to  yield 
data  of  a  similar  quality  to  data  obtained  using  a  custom-built  vehicle.  A  small  towed  platform 
was  designed  and  constructed  that  allowed  two  magnetometers,  two  EM61s,  and  a  GPS  to  be 
easily  slid  along  a  fiberglass  backbone  to  vary  the  vehicle-to-sensor  separation.  A  total  of  ten 
vehicles  were  tested,  ranging  in  size  from  a  small  ATV  up  through  a  Jeep  Wrangler,  but  the 
project’s  emphasis  was  on  side-by-side  UTVs  since  these  have  roll  cages,  shade,  ample  space  for 
electronics,  and  the  opportunity  for  weather-tight  enclosure,  whereas  small  ATVs  have  none  of 
these  things.  Six  UTVs  were  tested  using  the  adjustable  platform.  A  full  description  of  these 
vehicles  and  the  data  generated  is  contained  in  “MM-0605  Vehicle  Signature  Report  and 
Geophysical  Procedures  for  Vehicle  Signature  Measurement.doc”  submitted  to  the  Program 
Office  in  April  2007.  SAIC  selected  what  it  felt  was  the  best  compromise  of  low  signature,  cargo 
space,  terrain- handling  capability,  and  ride  -  the  aluminum- framed  Club  Car  XRT  1550  - 
purchased  it  with  SAIC  funds,  and  adapted  it  to  pull  SAIC’s  full-sized  VSEMS  mag/EM61 
towed  platform  (developed  under  ESTCP  project  MM-0208).  During  preliminary  testing  at  the 
small  Devens  test  site  in  Massachusetts,  we  found  that  this  vehicle  engenders  only  a  moderate 
signature  in  the  data,  and  that  this  signature  can  be  easily  removed  with  the  de-median  filter 
already  used  by  both  MTADS  and  VSEMS.  As  per  the  project  plan,  we  demonstrated  this 
vehicle,  towing  the  full-sized  VSEMS  platform,  at  APG,  to  measure  its  magnetic  signature  and 
evaluate  its  performance  on  a  near-real-world  site  with  a  greater  variety  of  terrain  than  the  small 
site  at  Devens. 

1.2  Objectives  of  the  Demonstration 

The  objectives  of  the  demonstration  were: 

•  To  acquire  survey  data  with  the  new  COTS  vehicle  over  a  site  that  has  realistic  terrain 
roughness,  and  compare  these  data  with  data  previously  acquired  over  the  same  site  using 
the  VSEMS  custom  aluminum- framed  buggy. 

•  To  apply  signature  removal  algorithms  to  the  data  and  judge  whether  there  is  a  significant 
difference  in  data  quality  attributable  to  the  COTS  vehicle  and  its  higher  bulk  magnetic 
signature. 
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•  To  evaluate  whether  there  are  significant  effects  in  signature-induced  artifacts  caused  by 
vehicle  pitch  and  roll  over  rough  terrain. 

•  To  evaluate  the  COTS  vehicle’s  terrain- handling  capability  and  ride  and  compare  them  to 
the  buggy’s. 

In  support  of  these  objectives,  we: 

•  Deployed  to  the  Standardized  UXO  Demonstration  Test  Site  at  Aberdeen  Proving 
Grounds,  MD 

•  Performed  octant  tests  (placed  the  system  in  a  clean  area  and  oriented  the  vehicle  and 
towed  platform  North,  Northeast,  East,  etc)  to  get  a  quick  snapshot  of  vehicle  signature  in 
several  configurations 

•  Surveyed  the  calibration  test  grid,  and  open  field  legacy  area  using  a  conservative 
configuration  (with  the  sensors  as  far  back  on  the  platform  as  they  can  possibly  go,  and 
without  the  steel  skidplate),  and  in  a  more  aggressive  configuration  (with  the  sensors  at 
their  nominal  15’  distance  from  the  vehicle,  and  with  the  skidplate) 

•  Examined  clean  areas  in  the  data  to  measure  bulk  vehicle  signature  and  perturbations  to 
that  bulk  signature  caused  by  vehicle  motion  on  uneven  terrain 

•  Varied  the  time  window  of  the  filtering  and  examined  the  effect  on  both  bulk  signature 
and  on  signal  and  amplitude  of  munitions  and  explosives  of  concern  (MEC)  targets 

1.3  Regulatory  Drivers 

The  primary  driver  is  the  continued  need  to  develop  tools  to  detect  MEC.  The  documented  use  of 
a  COTS  vehicle  to  tow  a  magnetometer  array  and  generate  high-quality  data  will  allow  more 
contractors  to  use  towed  magnetometer  arrays. 
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2  Technology 


2.1  Technology  Description 

2.1.1  Overview 

A  magnetometer’s  M/R3  omnidirectional  sensitivity  (where  M  is  magnetic  moment  in  amp- 
meters2,  and  R  is  the  distance  in  meters  between  detector  and  object)  makes  it  susceptible  to  the 
magnetic  signature  not  only  of  buried  MEC  but  also  of  a  towing  vehicle.  This  signature  changes 
primarily  with  the  vehicle’s  changing  orientation  relative  to  North,  and  changes  secondarily  as 
the  vehicle  pitches  and  rolls  about  the  horizontal  plane.  For  small  to  moderately- sized  vehicles, 
given  a  reasonable  distance  between  the  vehicle  and  the  towed  magnetometer  array,  the  vehicle’s 
signature  is  small  enough  and  changes  slowly  and  predictably  enough  with  orientation  relative  to 
North  that  the  signature  can  be  removed  with  the  same  filtering  techniques  used  to  remove  long- 
wavelength  geology,  diurnal  drift,  and  other  slowly-changing  effects.  As  vehicle  size  and  ferrous 
content  increases,  however,  small  fast  changes  in  vehicle  pitch  and  roll  begin  to  contribute 
substantially  to  vehicle  signature  in  a  way  that  defies  removal.  When  this  effect  begins  to 
dominate,  the  vehicle  is  too  large  and  /  or  the  sensors  are  too  close. 

2.1.2  Theory  of  Operation 

Any  ferrous  object  in  the  Earth’s  magnetic  field  has  a  magnetic  moment  induced  in  it.  The  field 
produced  by  that  moment  adds  vectorially  to  the  Earth’s  field.  The  magnitude  of  this  vector  sum 
is  measured  by  a  total  field  magnetometer.  During  data  processing,  the  magnitude  of  the  Earth’s 
field  is  then  subtracted  off,  leaving  the  anomalous  field.  From  sufficiently  far  away,  even  the 
complex  field  from  a  vehicle  resembles  a  point  dipole.  Grossly  speaking,  the  induced  dipole 
from  a  vehicle  will  tend  to  align  North,  and  will  tend  to  have  a  South-facing  positive  lobe,  and  a 
North-facing  negative  lobe.  This  means  that,  if  the  vehicle  towing  a  magnetometer  array  is 
headed  North,  the  magnetometers  will  be  in  the  positive  lobe  of  the  induced  dipole.  Conversely, 
if  the  vehicle  is  headed  South,  the  magnetometers  will  be  in  the  negative  lobe  of  the  anomaly 
created  by  the  induced  dipole.  This  is  depicted  in  the  conceptual  illustration  below.  Because  of 
the  directional  nature  of  the  vehicle’s  signature,  North-going  and  South-going  traverses  are 
affected  by  the  vehicle  differently,  and  contain  different  vehicle-induced  offsets.  Due  to  the  fact 
that  the  vehicle  is  an  extended  object,  not  a  point  dipole,  and  that  the  sensors  are  in  the  same 
plane  as  the  object  (not  above  the  object  as  with  detection  of  a  buried  object),  the  “in  the 
negative  lobe  of  the  anomaly”  explanation  is,  in  fact,  a  heuristic  simplification,  and  in  practice, 
different  vehicles  may  have  the  “North  positive”  and  “South  negative”  rules  flipped  around. 
However,  the  conceptual  illustration  is  extremely  useful  in  that  it  does  help  one  to  visualize  why 
the  measured  vehicle  signature  changes  in  passes  of  opposite  direction,  and  how  this  generates 
the  streaks  that  one  sees  in  vehicle-towed  magnetometer  data. 
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Mags  Behind  a 
Vehicle  Heading 
South  are  in  the 
Negative  Lobe 


Mags  Behind  a 
Vehicle  Heading  m 
North  are  in  the 
Positive  Lobe 


Figure  1:  Conceptual  illustration  depicting  how  survey  direction  affects  magnetic  signature 


The  image  below  on  the  left  shows  data  obtained  with  a  Yamaha  Rhino  (one  of  several  UTVs 
tested),  with  the  magnetometers  20’  behind  the  vehicle.  The  area  was  surveyed  in  a  racetrack- 
style  pattern,  with  the  traverses  on  the  left  side  of  the  image  obtained  in  North-going  passes,  and 
the  traverses  on  the  right  side  of  the  image  obtained  in  South-going  passes.  The  left  side  of  the 
image  (the  North-going  passes)  is  lighter  -  has  a  higher  average  value  -  than  the  right  side  of  the 
image  (the  South-going  passes).  For  this  vehicle,  at  this  20’  sensor-to-vehicle  separation,  the 
difference  between  North  and  South-going  traverses  is  about  40  nT.  This  value  is  calculated  by 
defining  small  areas  of  interest  in  each  of  the  North-going  and  the  South-going  sides,  each 
containing  no  other  magnetic  anomalies,  averaging  the  data  in  each  area,  and  then  subtracting  the 
two  averaged  values.  On  most  real  surveys,  a  racetrack  pattern  is  not  used,  and  instead,  adjacent 
traverses  are  usually  acquired  in  opposite  directions.  In  this  case,  we  divide  the  data  into  North¬ 
going  and  South-going  sets,  select  a  single  area  of  interest  that  contains  some  of  each  set, 
independently  accumulate  the  average  values  of  North-going  and  South-going  traverses,  and 
subtract  them.  In  either  case,  we  call  this  value  the  “bulk  signature.”  While  we  often  use  the 
terms  “bulk  signature”  and  “signature”  interchangeably,  the  “signature”  is  a  vector  quantity 
whose  magnitude  is  different  depending  on  the  vehicle’s  orientation  and  the  point  of 
measurement,  whereas  the  “bulk  signature”  is  a  single  number  representing,  for  a  fixed  vehicle- 
to-sensor  distance,  the  worst-case  difference  in  sensor  readings  taken  in  opposite  directions.  It  is 
usually  clear  from  the  context  which  is  meant. 


As  long  as  the  directionally-dependent  signature  is  the  only  dominant  effect,  it  changes  slowly 
and  can  be  easily  filtered  out.  In  the  images  below,  the  left  image  is  unfiltered,  and  the  right 
image  has  had  a  de-median  filter  with  a  six-second  time  window  applied  to  remove  the  vehicle 
signature  and  other  long-wavelength  effects.  Independent  analysis  by  both  SAIC  and  by  Dan 
Steinhurst  of  Nova  Research  (the  operator  of  vehicular  MTADS)  have  shown  that  this  sort  of  de¬ 
median  filter  does  a  good  job  at  removing  slowly-changing  vehicle  effects,  along  with  diurnal 
drift,  long-wavelength  geology,  and  small  inter-sensor  offsets  that  create  lines  along  the  direction 
of  travel  sometimes  informally  referred  to  as  “corduroy.” 
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Figure  2:  Sample  unfiltered  (left)  and  filtered  (right)  data  taken  with  one  of  several  UTVs  tested  (a  Yamaha 
Rhino,  mags  20’  back)  at  the  Devens  Test  Site,  clearly  showing  the  heading-dependent  effect  of  vehicle 
signature.  The  image  on  the  left  has  approximately  40nT  of  bulk  signature.  The  image  on  the  right  has  had 
bulk  signature  removed  through  de-median  filtering.  Image  scale  is  +-  50nT. 

The  key  here  is  the  phrase  “slowly-changing.”  Although  the  vehicle’s  signature  varies  primarily 
with  yaw  (orientation  relative  to  North),  it  also  changes  with  vehicle  roll  and  pitch.  As  terrain 
roughness  increases,  roll  and  pitch  increase,  creating  fast-changing  perturbations  in  the  data  that 
the  de-median  filter  can’t  take  out.  These  artifacts  show  up  as  a  “clouding”  of  the  data 
(essentially  a  disturbance  in  an  otherwise  smooth  background).  An  example  of  “cloudy”  data  is 
shown  in  the  figure  on  the  right,  below.  Note  that  both  of  these  images  have  been  de-median 
filtered  using  a  six-second  window,  the  effect  of  which  has  been  to  remove  the  azimuthally- 
varying  vehicle  signature  (that  is,  there  is  no  obvious  light/dark  bias  in  either  image  resulting 
from  the  difference  in  signature  North  to  South),  but  the  image  on  the  right  clearly  has  artifacts 
that  the  image  on  the  left  does  not.  These  are  due  to  pitch  and  roll-based  changes  in  vehicle 
signature  caused  by  localized  terrain  roughness.  Short  of  putting  the  vehicle  on  a  jig  that  allows 
for  careful  roll  and  pitch  of  the  entire  vehicle  and  taking  a  full  set  of  measurements,  it  is  difficult 
to  quantify  and  predict  the  effect  of  pitch  and  roll.  Some  sites  are  simply  rougher  than  others,  and 
for  a  given  vehicle  and  sensor  distance,  some  sites  may  engender  more  of  these  sort  of  artifacts 
than  other  sites.  For  example,  if  one  had  to  survey  a  very  smooth  section  of  beach  on  an  island, 
and  it  was  prohibitively  expensive  to  transport  a  vehicle  onto  the  island,  and  if  a  steel  UTV  such 
as  a  Gator  was  available  on-island,  use  of  such  a  vehicle  might  be  plausible.  But  if  these  artifacts 
appear  in  the  data  at  a  magnitude  that  begins  to  interfere  with  detection  and  analysis  of  the 
targets  of  interest,  the  vehicle  signature  is  too  high  and  the  magnetometers  are  too  close. 
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Figure  3:  Data  taken  with  steel  John  Deere  Gator  tow  vehicle  showing  “cloudy”  artifacts  due  to  pitch  and  roll 
over  rough  terrain.  Left  image  is  with  sensors  20  feet  back,  right  image  is  with  sensors  10  feet  back.  Both 
images  are  de-median  filtered  and  shown  at  +-  50  nT.  The  terrain  at  the  northern  and  southwestern  edge  of 
the  Devens  site  is  rough,  causing  vehicle  pitch  and  roll  that  changes  the  magnetic  signature  seen  at  the 
sensors.  The  de-median  filter  takes  out  the  bulk  signature  at  both  sensor  separations,  but  with  the  sensors  at 
10’  the  effects  from  pitch  and  roll  become  substantial  and  are  not  taken  out  by  the  filter. 


2.1.3  Key  Design  Criteria 

The  key  design  criteria  are  the  trade-offs  between  the  vehicle’s  ferrous  content,  the  separation 
between  the  vehicle  and  the  magnetometers,  the  maneuverability  of  the  vehicle,  the  cargo 
carrying  capability  of  the  vehicle,  and  the  median  filter  window  length. 

2.1.4  Schematics,  Figures,  and  Layout 

The  figure  below  shows  the  Club  Car  1550XRT  towing  the  VSEMS  platform. 
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Figure  4:  Club  Car  1550XRT  towing  the  VSEMS  platform 


2.2  Technology  Development 

•  Fall  2006.  A  simple  two-magnetometer  and  two-EM61  towed  “TSEMS”  platform  was 
developed  to  allow  vehicle  signatures  to  be  evaluated  with  adjustable  vehicle-to-sensor 
separation.  Ten  vehicles  ranging  from  an  ATV  to  a  Jeep  were  procured  through 
borrowing  or  renting,  and  tested  at  the  Devens  MA  test  site  over  a  four-month  period. 

•  Spring  2007.  Data  analyzed.  Aluminum-framed  vehicles  (Club  Car  and  Bobcat)  were 
determined  to  have  best  signature  characteristics  of  all  the  UTVs.  Verified  that  a  de¬ 
median  filter  did  a  good  job  removing  the  already  low  vehicle  signature.  Unfortunately 
these  particular  aluminum-framed  vehicles  do  not  have  fully-independent  rear 
suspension.  This  resulted  in  a  very  stiff  ride.  We  had  concerns  that  the  ride  was  so  stiff 
that  it  might  make  the  vehicle  a  poor  choice  for  day-in  day-out  DGM,  but  the 
manufacturer  told  us  that  the  suspension  was  designed  to  have  6001bs  of  load  in  the  dump 
bed.  We  felt  that,  once  the  vehicle  was  equipped  with  all  the  electronics  (which  includes 
seven  car  batteries),  the  suspension  would  be  sufficiently  pre-loaded. 

•  Fall  2007.  Procured  demonstrator  model  of  2007  Club  Car  for  testing  (the  previous  one 
tested  was  2004  model)  to  be  absolutely  certain  that  the  vehicle  we  were  about  to 
purchase  would  have  signature  similar  to  the  one  tested. 

•  Winter  2008.  Purchased  2007  Club  Car  1550XRT,  with  enclosure  and  heat,  using  SAIC 
funds. 

•  Spring  2008.  Performed  final  adaptation  of  Club  Car  1550XRT  to  pull  VSEMS  trailer. 
This  involved  removing  the  dump  bed,  mounting  all  VSEMS  electronics  in  a  clamshell 
on  the  back  deck,  outfitting  with  alloy  wheels,  demagnetizing  tire  beads,  and  adapting  the 
trailer  hitch  for  the  disparate  vehicle  /  platform  height.  Performed  preliminary  testing  of 
vehicle  at  Devens  test  site.  The  pre-loading  of  the  suspension  from  the  weight  of  the 
seven  car  batteries  made  the  ride  quality  of  the  fully-loaded  vehicle  acceptable. 

2.3  Advantages  and  Limitations  of  the  Technology 

Towed  array  magnetometry  is  most  directly  applicable  to  large,  vehicular ly- hospitable  sites  that 
contain  major  caliber  air-dropped  munitions  (e.g.,  2501b  bombs)  and  do  not  have  geology  with 
iron-bearing  soil.  If  munitions  of  interest  are  between  60mm  and  155mm,  published  reports 
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indicate  that  magnetometers  and  EM61s  have  a  similar  performance  envelope.  For  munitions 
smaller  than  60mm  with  little  ferrous  content,  published  reports  indicate  that  the  performance 
advantage  tips  toward  the  EM61.  As  terrain  becomes  more  rugged  and  less  hospitable  to 
vehicular  surveys,  vehicle-towed  arrays  themselves  become  less  applicable,  regardless  of  the 
choice  of  sensor. 

The  alternative  technology  to  a  COTS  vehicle  is  development  of  a  custom  low-ferrous  vehicle 
with  an  aluminum  frame.  Because  of  the  high  (>$100k)  cost,  there  are  only  three  of  these 
vehicles  in  existence,  all  built  by  Chenowth  Racing.  These  vehicles,  sometimes  incorrectly 
referred  to  as  “dune  buggies”  (they  are  not;  they  are  off-road  racing  vehicles)  have  worked  well 
(the  MTADS  vehicle  is  still  in  use,  and  the  STOLS  /  VSEMS  vehicle  was  retired  only  recently), 
but  unlike  a  COTS  UTV,  they  do  not  have  four-wheel  drive.  They  are  difficult  to  make 
weatherproof  (unlike  a  COTS  vehicle  with  an  enclosure).  Further,  even  the  custom  $100k 
aluminum-framed  dune  buggies  have  ferrous  components  (the  internal  engine  and  transaxle 
components,  for  example,  are  steel),  and  because,  for  practical  reasons,  you  can’t  tow  sensors 
infinitely  far  back,  even  these  vehicles  have  a  non-zero  bulk  magnetic  signature  at  practical 
towing  distances  of  15  or  20  feet. 

Although  signature  is  important,  the  vehicle’s  capacity  to  carry  needed  electronics  and  the 
vehicle’s  terrain  handling  ability  are  also  critical.  Thus,  this  project  involved  testing  a  total  of  ten 
new  vehicles  and  finding  the  best  tradeoff  between  vehicle  size,  vehicle  signature,  terrain 
handling  capacity,  and  vehicle-to-sensor  separation.  These  results  are  contained  in  “MM-0605 
Vehicle  Signature  Report  and  Geophysical  Procedures  for  Vehicle  Signature  Measurement.doc” 
previously  submitted  to  the  ESTCP  Program  Office.  The  primary  conclusions  of  the  research 
were: 

•  Due  to  their  small  size,  ATVs  don’t  have  much  ferrous  metal  and  thus  have  low  magnetic 
signatures,  but  their  lack  of  roll  cage,  inability  to  be  made  weather-tight,  small  storage 
space,  and  straddle- style  seating  make  them  less  than  optimal  tow  vehicles  for  real-world 
production  DGM. 

•  In  contrast,  side-by-side  UTVs  offer  a  roll  cage,  ample  storage  space,  upright  seating,  and 
the  possibility  of  a  weather-enclosed  cabin.  However,  the  resulting  increase  in  ferrous 
metal  creates  a  non-trivial  magnetic  signature. 

•  Of  the  UTVs  tested,  the  vehicle  with  the  smallest  signature  was  an  aluminum- framed 
vehicle  jointly  designed  by  Bobcat,  Club  Car,  and  Husqvarna.  Because  of  its  aluminum 
frame  and  consequently  reduced  ferrous  mass,  this  vehicle  had  a  substantially  smaller 
signature  than  any  other  full-sized  UTV. 

•  We  tested  a  Club  Car  1550XRT  aluminum- framed  vehicle  with  its  steel  dump  bed 
removed,  its  steel  wheels  replaced  with  alloy  wheels,  and  its  tire  beads  demagnetized,  and 
found  it  to  be  a  very  capable  survey  vehicle  with  a  very  low  magnetic  signature  -  not 
quite  as  low  as  the  Chenowth  buggy,  but  low  enough  so  that  the  remaining  vehicle 
signature  could  be  easily  removed  with  a  de-median  filter. 

•  Because  SAIC’s  VSEMS  system  concurrently  collects  both  EM61  data  and 
magnetometer  data,  we  selected  a  Club  Car  1550XRT  with  a  diesel  engine  instead  of  a 
gas  engine.  The  diesel  engine  has  a  slightly  higher  magnetic  signature  than  the  gas  engine 
due  to  its  ferrous  block,  but  the  diesel  engine  has  a  lower  RF  noise  output  than  the  gas 
engine  due  to  the  fact  that  diesel  engines  lack  an  ignition  system.  The  version  of  this 
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vehicle  equipped  with  a  gas  engine  had  a  lower  bulk  magnetic  signature,  but  more  noise 
on  the  EM61s.  Note  that  SAIC  purchased  its  diesel  1550XRT  as  capital  equipment;  the 
purchase  of  this  vehicle  was  not  charged  to  project  MM-0605. 
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3  Performance  Objectives 

The  objective  of  the  project  as  a  whole  was  to  test  a  variety  of  COTS  vehicles,  select  one  with  a 
low  magnetic  signature,  configure  it  as  a  tow  vehicle,  use  it  on  a  real  towed  magnetometer 
survey,  and  demonstrate  that,  with  filtering,  the  resulting  data  are  comparable  to  data  acquired 
using  the  $100k  custom  Chenowth-built  aluminum  vehicle  employed  by  VSEMS  and  similar  to 
that  employed  by  MTADS.  Thus  the  objectives  of  the  demonstration  were: 

•  To  acquire  survey  data  with  the  new  COTS  vehicle  over  a  site  that  has  realistic  terrain 
roughness,  and  compare  these  data  with  data  previously  acquired  over  the  same  site  using 
the  VSEMS  custom  aluminum- framed  buggy. 

•  To  apply  signature  removal  algorithms  to  the  data  and  judge  whether  there  is  a  significant 
difference  in  data  quality  attributable  to  the  COTS  vehicle  and  its  higher  bulk  magnetic 
signature. 

•  To  evaluate  whether  there  are  significant  effects  in  signature-induced  artifacts  caused  by 
vehicle  pitch  and  roll  over  rough  terrain. 

•  To  evaluate  the  COTS  vehicle’s  terrain- handling  capability  and  ride  and  compare  them  to 
the  buggy’s. 

In  support  of  these  objectives,  we: 

•  Deployed  to  the  Standardized  UXO  Demonstration  Test  Site  at  Aberdeen  Proving 
Grounds,  MD 

•  Performed  octant  tests  (placed  the  system  in  a  clean  area  and  oriented  the  vehicle  and 
towed  platform  North,  Northeast,  East,  etc)  to  get  a  quick  snapshot  of  vehicle  signature  in 
several  configurations 

•  Surveyed  the  calibration  test  grid,  and  open  field  legacy  area  using  a  conservative 
configuration  (with  the  sensors  17’  behind  the  vehicle  -  as  far  back  on  the  platform  as 
they  can  possibly  go  -  and  without  the  steel  skidplate),  and  in  a  more  aggressive 
configuration  (with  the  sensors  at  their  nominal  15’  distance  from  the  vehicle,  and  with 
the  skidplate) 

•  Examined  clean  areas  in  the  data  to  measure  bulk  vehicle  signature  and  perturbations  to 
that  bulk  signature  caused  by  vehicle  motion  on  uneven  terrain 

•  Varied  the  time  window  of  the  filtering  and  examine  the  effect  on  both  bulk  signature  and 
on  signal  and  amplitude  of  MEC  targets 

The  performance  objectives  are  listed  in  the  table  below. 
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Table  1:  Performance  Objectives 


Performance 

Objective 

Metric 

Data  Required 

Success  Criteria 

Quantitative  Performance  Objectives 

Small  to  moderate 
directional  effect 
of  vehicle 
magnetic  signature 
in  unfiltered  data 

Bulk  magnetic  signature 
in  unfiltered  data 

•  Static  octant  data 

•  Dynamic  survey  data 

<  25  nT 

Small  directional 
effect  of  vehicle 
magnetic  signature 
in  filtered  data 

Bulk  magnetic  signature 
in  filtered  data 

•  Dynamic  survey  data 

<2nT 

Small  effect  of 
filtering  on 
amplitude  and  size 
of  targets 

Amplitude  above 
background  and  size 
(full  width  at  half 
maximum  value) 

•  Dynamic  survey  data 

•  nine  targets 

<15%  difference 
(filtered  vs  unfiltered) 
in  amplitude  and  size 

Small  pitch  and 
roll  effect  of 
vehicle  magnetic 
signature 

Magnitude  of  artifacts  in 
survey  data  over  smooth 
and  rough  terrain 

•  GPS  locations  of 
smooth  and  rough 
areas  within  APG 
open  field  legacy 
area 

•  Dynamic  survey  data 
over  smooth  terrain 

•  Dynamic  survey  data 
over  rough  terrain 

<  2  nT  over  smooth 

areas 

<  5  nT  over  rough 
areas 

Qualitative  Performance  Objectives 

Vehicle  Ride 

Operator  Observations 

•  Having  operator 
drive  the  vehicle 
during  a 

representative  survey 

Ride  is  acceptable  - 
not  so  stiff  that 
operator  is  beaten  up 

Vehicle  Terrain- 
Handling 

Capability 

Operator  Observations 

•  Having  operator 
drive  the  vehicle 
during  a 

representative  survey 

Vehicle  handles  terrain 
at  least  as  well  as 
Chenowth  vehicle 

All  of  the  above  objectives  were  met.  The  demonstration  showed  that  a  COTS  vehicle  can  collect 
data  that,  when  filtered,  are  comparable  with  data  acquired  with  a  custom-built  vehicle. 
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4  Site  Description 

4.1  Site  Selection 

The  demonstration  was  conducted  at  the  Standardized  UXO  Demonstration  Test  Site  in 
Aberdeen,  MD.  This  site  was  selected  because  we  had  surveyed  it  before  with  VSEMS  towed  by 
a  Chenowth  vehicle,  and  because  the  site  has  sections  that  are  very  rough,  enabling  direct 
examination  of  the  degree  to  with  a  rough  site  engenders  pitch  and  roll  changes  in  the  vehicle 
which  in  turn  creates  perturbations  in  the  signature  that  creates  artifacts  in  the  filtered  data. 

4.2  Site  History 

The  APG  and  YPG  sites  were  established  in  1999  to  provide  a  standard  demonstration  area  for 
emerging  MEC  detection-related  technologies.  The  Standardized  Technology  Demonstration 
Test  Site  in  Aberdeen  MD  has  a  calibration  grid,  a  blind  test  grid,  and  an  open  field.  The  open 
field  itself  contains  a  direct  fire,  indirect  fire,  and  legacy  area;  it  is  the  legacy  area  that  we  will 
concentrate  on.  The  calibration  lane  and  blind  test  grid  have  no  surface  features  of  concern.  The 
open  field  site  is  generally  flat  with  a  low  area  that  is  wet  during  a  portion  of  the  year,  a  tree 
stump  area,  and  a  section  of  gravel  road. 

4.3  Site  Geology 

The  geology  at  the  APG  site  poses  no  challenges  to  the  magnetometers. 

4.4  Munitions  Contamination 

MEC  at  the  site  ranges  from  20mm  projectiles,  40mm  grenades,  submunitions,  60  -  81mm 
mortars,  2.75  inch  rockets,  105-155mm  projectiles,  and  2501b  bombs. 

5  Test  Design 

5.1  Conceptual  Experimental  Design 

As  described  above,  during  the  project,  we  tested  ten  vehicles  ranging  in  size  from  a  small  ATV 
to  a  Jeep  as  candidates  to  replace  the  custom  Chenowth  tow  vehicle.  We  concentrated  on  side- 
by-side  UTVs  as,  unlike  ATVs,  these  have  roll  cages,  shade,  upright  seating,  ample  cargo  space, 
and  the  ability  to  be  enclosed  against  weather  (ATVs  have  none  of  these).  We  selected  the  Club 
Car  1550XRT  as  the  vehicle  we  felt  to  be  the  best  compromise  of  vehicle  signature,  terrain 
handling,  cargo  capacity,  and  ride  quality.  Tests  were  conducted  on  a  small  area  in  Devens,  MA 
with  limited  terrain  variation.  The  goal  of  the  demonstration  was  to  test  the  vehicle  towing  the 
full-sized  VSEMS  platform  on  a  site  that  was  larger  than  the  Devens  site  and  had  greater 
variation  of  terrain  then  the  Devens  site,  and  to  evaluate  1)  what  the  vehicle’s  bulk  signature  is; 
2)  whether  that  signature  can  be  effectively  removed  through  filtering  without  adversely 
affecting  signals  from  buried  MEC;  3)  how  the  unfiltered  and  filtered  data  compare  with  data 
taken  with  the  custom  Chenowth  vehicle  in  2006;  4)  whether  the  signature  is  large  enough  that 
there  are  significant  pitch  and  roll  effects  that  can’t  be  removed  with  filtering;  and  5)  whether  the 
vehicle  is  usable  for  real-world  DGM  in  terms  of  terrain-handling  capability  and  ride  quality. 
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5.2  Site  Preparation 

The  APG  site  needs  no  preparation.  We  did,  however,  acquire  GPS-based  geodetic  coordinates 
to  outline  rough  and  smooth  areas  so  that  we  can  later  correlate  pitch  and  roll-based  artifacts  in 
the  data  with  varying  site  topography.  We  did  this  post-survey  since,  at  that  point,  the  driver 
knew  where  the  smooth  and  rough  areas  were. 

5.3  System  Specification 

We  deployed  the  selected  COTS  vehicle  -  a  Club  Car  1550XRT  -  adapted  to  tow  the  VSEMS 
towed  platform.  The  VSEMS  towed  platform  is  purpose-built  to  host  both  magnetometers  and 
EM61  coils  in  a  low-noise  environment.  It  is  constructed  of  carbon  fiber  and  utilizes  springs  and 
air  bags  to  absorb  and  damp  out  terrain- induced  sensor  motion.  The  platform’s  tires  have  had 
their  steel  beads  removed,  and  the  platform’s  wheels  contain  no  rotating  ferrous  metal.  At  the 
back  of  the  platform  are  five  Geometries  822A  cesium  vapor  total  field  magnetometers,  1’  (0.30 
meters)  above  ground,  spaced  0.5  meters  apart.  The  magnetometer  data  acquisition  hardware 
utilizes  the  interleaving  design  developed  under  ESTCP  project  MM-0208  that  acquires 
magnetometer  data  at  75  Hz  between  EM61  pulses.  Five  1  x  Vi  meter  coils  are  employed,  with 
the  short  axis  oriented  cross-track.  EM61  output  occurs  at  10  Hz.  Although  EM61  data  is  not 
required  for  this  project,  since  the  capability  to  concurrently  collect  it  is  integral  with  VSEMS, 
we  collected  it.  Three  Trimble  MS750  RTK  GPS  receivers  were  employed,  with  one  antenna 
over  the  center  magnetometer,  one  over  the  center  of  the  EM61  array,  and  one  on  the  roof  of  the 
vehicle  providing  input  to  a  track  guidance  system. 

Adaption  of  the  Club  Car  vehicle  involved  transferring  all  sensor,  positioning,  data  logging,  and 
power  electronics  from  the  VSEMS  Chenowth  vehicle.  This  was  accomplished  by  removing  the 
Club  Car’s  steel  dump  bed  (conveniently,  the  largest  piece  of  ferrous  metal  on  the  vehicle), 
replacing  it  with  a  COTS  plastic  rooftop  car  carrier,  and  mounting  the  power,  electronics,  and 
computing  systems  inside.  A  computer  monitor  and  power  control  switches  were  mounted  inside 
the  vehicle. 

Note  that  part  of  the  project  philosophy  was  that  the  COTS  vehicle  remain  as  COTS  as  possible. 
That  is,  easily  unbolt-able  steel  components  such  as  the  steel  dump  bed  were  removed,  and  the 
steel  wheels  were  replaced  with  the  alloy  wheels  directly  available  from  the  manufacturer,  but 
we  did  not  have  any  custom  non-ferrous  parts  fabricated.  Thus  the  roll  cage,  the  rear  subframe 
and  rear  trailing  arms,  and  the  trailer  hitch  are  all  steel.  In  addition,  we  ordered  the  two  skidplates 
that  are  available  from  the  manufacturer  to  protect  the  vehicle’s  underside  -  particularly  the  front 
and  rear  driveshafts  -  from  rocks.  The  front  skidplate  is  aluminum,  but  the  rear  skidplate  is  steel, 
and  is  not  simply  a  flat  steel  plate;  it  is  a  contoured  part  made  up  of  many  welded  steel  sections. 
It  certainly  is  possible  to  fabricate  the  same  design  out  of  aluminum,  but  it  stretches  the  project’s 
COTS  philosophy.  As  such,  we  procured  the  steel  rear  skidplate,  and  tested  the  vehicle  both  with 
and  without  it. 

Use  of  VSEMS  also  includes  a  reference  magnetometer  employed  to  measure  diurnal  variations 
in  the  Earth’s  field.  We  have  not  made  frequent  use  of  the  reference  magnetometer  for  several 
years,  instead  relying  on  use  of  a  de-median  filter  to  remove  all  long-wavelength  effects 
including  diurnal  drift,  geology,  and  vehicle  signature.  However,  for  this  demonstration,  we 
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wished  to  be  able  to  separate  these  effects,  so  we  deployed  the  reference  magnetometer  and 
downloaded  its  data  at  the  end  of  each  day. 

5.3.1  Operating  Parameters  for  the  Technology 

As  said  above,  the  largest  “operating  parameter”  was  the  selection  of  the  vehicle  (our  Vehicle 
Signature  Report  contains  extensive  detail  on  the  signatures  of  all  ten  vehicles  we  tested,  each  at 
four  different  sensor  separations).  We  went  to  APG  to  demonstrate  with  the  selected  vehicle.  We 
did,  however,  vary  the  vehicle-to-sensor  separation  by  moving  the  magnetometers  back  from 
their  nominal  15’  location  to  a  17’  location  at  the  furthest  point  on  the  back  of  the  survey 
platform.  In  addition,  we  varied  the  ferrous  mass  of  the  vehicle  by  attaching  and  removing  the 
COTS  ferrous  welded  steel  skidplate  intended  to  protect  the  driveshaft  and  other  undercarriage 
components.  As  per  the  demonstration  plan,  we  collected  octant  data  in  four  configurations  (15” 
with  skidplate  on;  15’  with  skidplate  off;  17’  with  skidplate  on;  17’  with  skidplate  off),  and 
collected  survey  data  in  the  two  endpoint  configurations  (15’  with  skidplate  on;  17’  with 
skidplate  off). 

5.4  Calibration  Activities 

There  were  no  calibration  activities  per  se.  However,  we  acquired  several  baseline  octant  test 
data  sets.  These  octant  test  data  sets  were  not  used  to  generate  any  metrics  used  in  the  evaluation; 
however  they  are  a  quick  and  effective  way  of  providing  a  snapshot  of  how  the  vehicle  signature 
changes  with  increasing  sensor  separation  and  with  the  addition  of  ferrous  metal  to  the  vehicle. 
We  positioned  the  magnetometer  sensor  array  on  a  magnetically  clean  area  near  the  test  pit  and 
oriented  the  vehicle  along  magnetic  North,  South,  East,  West,  and  the  four  bisected  directions, 
and  collected  30  seconds  of  data  in  each  orientation.  We  did  this  for  four  configurations  noted 
above. 

5.5  Data  Collection 

5.5.1  Scale 

The  octant  tests  were  performed.  A  5.5  acre  portion  of  the  legacy  area  of  the  open  field  was 
surveyed.  We  also  surveyed  the  calibration  grid. 

5.5.2  Sample  Density 

The  interleaving  hardware  in  VSEMS  outputs  magnetometer  data  at  75  Hz.  At  a  nominal  survey 
speed  of  two  meters  per  second  (approximately  4.6  mph),  this  creates  a  magnetometer  data  value 
every  2.7  cm  along  the  ground.  The  GPS  outputs  at  10  Hz.  Adjacent  survey  lines  were  spaced  to 
overlap  the  outer  magnetometer.  These  parameters  were  also  used  to  acquire  VSEMS  data  at 
APG  in  2006 

5.5.3  Quality  Checks 

Equipment  was  warmed  up  for  five  minutes  prior  to  survey  operations.  In-vehicle  software 
displays  both  numeric  and  graphic  representations  of  magnetometer,  GPS,  and  EM61  output. 
The  RTK  GPS  outputs  a  fix  quality  of  3  for  a  centimeter-level  fixed  integer  solution.  When  the 
RTK  solution  is  3,  the  survey  screen’s  background  is  gray.  For  an  RTK  solution  of  2  (“RTK 
float”),  the  survey  screen  is  colored  yellow,  indicating  to  the  operator  that  RTK  is  in  the  process 
of  re-initializing.  The  operator  knows  that  if  the  screen  does  not  go  gray  (indicating  RTK  3)  very 
quickly,  he  should  stop  driving  the  vehicle.  Fix  qualities  of  1  (autonomous  GPS,  indicating  that 
the  rover  GPS  in  the  vehicle  is  not  receiving  corrections  from  the  base  station),  4  (differential 
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non-RTK  GPS),  and  0  (not  enough  satellites  to  generate  a  fix)  result  in  the  screen  being  colored 
bright  red,  telling  the  operator  to  stop  driving  immediately  since  useless  data  are  being  acquired. 
Non-RTK  data  was  only  seen  when  driving  along  the  western  wood  line.  The  updates  from  all 
five  magnetometers  are  displayed  waterfall-style,  enabling  the  operator  to  continuously  monitor 
their  performance  while  driving. 

5.5.4  Data  Summary 

We  collected  static  octant  data: 

•  With  the  magnetometers  at  their  current  15’  distance  from  the  vehicle  and  with  the 
vehicle’s  steel  skidplate  attached 

•  With  the  magnetometers  at  their  current  15’  distance  from  the  vehicle  and  with  the 
vehicle’s  steel  skidplate  removed 

•  With  the  magnetometers  moved  back  an  additional  two  feet  and  with  the  vehicle’s  steel 
skidplate  attached 

•  With  the  magnetometers  moved  back  an  additional  two  feet  and  with  the  vehicle’s  steel 
skidplate  removed 

We  collected  dynamic  survey  data  over  a  5.5  acre  section  of  the  legacy  area  and  over  the 
calibration  grid: 

•  With  the  magnetometers  at  their  current  15’  distance  from  the  vehicle  and  with  the 
vehicle’s  steel  skidplate  attached 

•  With  the  magnetometers  moved  back  an  additional  two  feet  and  with  the  vehicle’s  steel 
skidplate  removed 

These  data  reside  at  SAIC  in  Waltham,  MA,  on  the  server,  in  an  ASCII  comma-delimited  format 
(easting,  northing,  sensor_number,  line_number,  time,  nT)  format. 

5.6  Validation 

No  digging  was  performed  on  this  project. 
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6  Data  Analysis  and  Products 

The  overall  plan  was  to: 

•  Correct,  pre-process,  process,  and  image  the  data  in  the  standard  fashion  for  towed  array 
magnetometer  data  (with  a  de- median  filter  to  remove  the  directionally-changing  vehicle 
signature) 

•  Evaluate  the  filtered  and  the  unfiltered  data  by  dividing  the  data  into  sets  of  directionally- 
like  data  (e.g.,  North-going  and  South-going)  and  averaging  data  in  areas  of  interest  that 
contain  no  anomalies 

•  Vary  the  time  window  of  the  de-median  filter  to  determine  the  value  that  takes  out 
vehicle  signature  without  changing  size  or  amplitude  of  buried  MEC  by  more  than  15% 

•  Compare  the  bulk  signatures  of  the  unfiltered  and  filtered  data  with  those  obtained  with 
the  Chenowth  vehicle  in  2006 

•  Evaluate  the  presence  and  strength  of  artifacts  caused  by  perturbations  in  the  signature 
due  to  changes  in  vehicle  pitch  and  roll 

6.1  Preprocessing 

Data  correction  was  performed  on  both  the  GPS  data  and  the  magnetometer  data.  All  GPS  data 
were  displayed  graphically  and  examined  by  a  trained  analyst.  Any  GPS  readings  that  was  not  of 
fix  quality  3  (centimeter-level  RTK  fixed  integer  solution)  was  flagged  on  the  screen.  Individual 
GPS  jumps  were  corrected  via  interpolation.  Any  large  sections  of  non-RTK3  data  were  thrown 
out.  The  sensor  array’s  orientation  relative  to  North  was  then  calculated  by  first  temporarily 
smoothing  the  GPS  values,  and  then  using  adjacent  GPS  readings  to  calculate  platform  heading. 
The  sensor  array’s  orientation  was  then  calculated  as  normal  to  the  platform’s  heading.  The 
sensor  data  were  also  viewed  by  an  analyst  to  flag  any  noisy  data  that  made  it  passed  the  real¬ 
time  quality  check  in  the  vehicle.  Spurious  spikes  were  removed  via  a  median  filter. 

Next,  a  notch  filter  was  run  on  the  magnetometer  data.  Because  VSEMS  acquires  concurrent 
mag  and  EM61  data,  the  magnetometer  sampling  occurs  at  the  EM61’s  75  Hz  pulse  repetition 
rate.  At  75  Hz,  the  ubiquitous  60  Hz  hum  from  ambient  electrical  activity  aliases  flawlessly  at  15 
Hz.  A  de-spiking  median  filter  was  first  applied  to  the  time-series  magnetometer  data  on  each 
line  to  remove  spurious  values.  Then,  a  notch  filter  was  applied  to  the  magnetometer  data  to 
remove  the  15  Hz  aliased  signal. 

Data  from  the  reference  magnetometer  was  read  in,  time-correlated  with  the  survey  data,  and 
subtracted  to  remove  diurnal  drift.  This  drift  is  also  corrected  by  the  background  leveling  step 
below,  but  we  wish  to  correct  for  it  individually  because  we  will  need  unleveled  data  to  calculate 
the  bulk  signature  before  removal  by  the  de- median  filter. 

A  de-median  filter  with  a  six-second  window  is  nominally  applied  to  the  magnetometer  data, 
separately  for  each  magnetometer,  to  determine  a  background  value.  This  value  is  then 
subtracted  off  from  the  data,  resulting  in  dynamic  background  leveling.  This  also  removes 
geology  and  both  remnant  and  induced  vehicle  magnetic  signature.  The  length  of  the  median 
filter  window  was  varied  in  a  separate  step. 
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The  magnetometer  data  and  the  geolocation  data  are  then  combined.  VSEMS  utilizes  the  GPS’ 
1PPS  output  to  trigger  the  acquisition  of  one  second’s  worth  of  magnetometer  data.  In  this  way, 
the  whole-second-aligned  GPS  solution  corresponding  to  the  1PPS  is  guaranteed  to  correspond 
with  the  first  set  of  75  readings  in  the  one-second  data  block.  This  creates  perfectly  synchronized 
magnetometer/GPS  data  that  does  not  need  to  be  lag-corrected.  The  array  orientation  calculated 
from  the  GPS-synthesized  platform  heading  is  then  combined  with  the  array  geometry  to 
calculate  the  location  of  each  of  the  five  magnetometers  for  each  of  the  75  updates  within  each 
second. 

The  above  preprocessing  occurs  in  our  own  software.  Data  are  then  written  out  in  an  ASCII 
Geosoft  Oasis-importable  format. 

6.2  Target  Selection  for  Detection 

Target  selection  for  detection  was  not  performed.  As  per  the  demonstration  plan,  we  looked  at 
bulk  signature  and  at  the  effects  of  filtering. 

6.3  Parameter  Estimates 

Not  applicable. 

6.4  Classifier  and  Training 

Not  applicable. 

6.5  Data  Products 

ASCII  comma-delimited  files  as  described  in  section  5.5.4  were  produced.  These  files  were 
imported  into  Geosoft  Oasis  Montaj,  and  gridded  data  and  maps  were  produced. 
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7  Performance  Assessment 


7.1  Performance  Criteria 

Performance  criteria  from  the  Demonstration  Plan  are  listed  the  table  below. 


Performance 

Objective 

Metric 

Data  Required 

Success  Criteria 

Criteria 

Met? 

Quantitative  Performance  Objectives 

Small  to  moderate 
directional  effect 
of  vehicle 
magnetic  signature 
in  unfiltered  data 

Bulk  magnetic 
signature  in 
unfiltered  data 

•  Static  octant  data 

•  Dynamic  survey 
data 

<  25  nT 

Yes 

Small  directional 
effect  of  vehicle 
magnetic  signature 
in  filtered  data 

Bulk  magnetic 
signature  in 
filtered  data 

•  Dynamic  survey 
data 

<2nT 

Yes 

Small  effect  of 
filtering  on 
amplitude  and  size 
of  targets 

Amplitude 

above 

background 
and  size  (full 
width  at  half 
maximum 
value) 

•  Dynamic  survey 
data 

•  Nine  targets 

<  15% 
difference 
(filtered  vs 
unfiltered)  in 
amplitude  and 
size 

Yes 

Small  pitch  and 
roll  effect  of 
vehicle  magnetic 
signature 

Magnitude  of 
artifacts  in 
survey  data 
over  smooth 
and  rough 
terrain 

•  GPS  locations  of 
smooth  and  rough 
areas  within  APG 
open  field  legacy 
area 

•  Dynamic  survey 
data  over  smooth 
terrain 

•  Dynamic  survey 
data  over  rough 
terrain 

<  2  nT  over 
smooth  areas 

<  5  nT  over 
rough  areas 

Yes 

Qualitative  Performance  Objectives 

Vehicle  Ride 

Operator 

Observations 

•  Having  operator 
drive  the  vehicle 
during  a 
representative 
survey 

Ride  is 

acceptable  -  not 
so  stiff  that 
operator  is 
beaten  up 

Yes 

Vehicle  Terrain- 
Handling 

Capability 

Operator 

Observations 

•  Having  operator 
drive  the  vehicle 
during  a 
representative 
survey 

Vehicle  handles 
terrain  at  least 
as  well  as 
Chenowth 
vehicle 

Yes 
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7.2  Performance  Confirmation  Methods 

We  used  the  COTS  Club  Car  1550XRT  vehicle  to  tow  the  VSEMS  magnetometer  array.  The 
bulk  signature  was  measured  by  surveying  a  5.5  acre  section  of  the  Legacy  Area  in  two  different 
configurations.  It  was  also  measured  by  performing  octant  tests. 

The  effect  of  filtering  was  measured  by  performing  a  de-median  filter  in  Oasis,  varying  the  time 
window  length  of  the  filter,  and  comparing  the  amplitude  of  nine  objects. 

The  effect  of  vehicle  pitch  and  roll  in  the  data  was  measured  by  using  GPS  landmark  data  to 
identify  rough  areas  in  the  survey  site,  selecting  an  area  of  interest  in  a  rough  section,  and 
performing  a  statistical  comparison  with  an  area  of  interest  in  a  smooth  section. 

The  vehicle’s  handling  and  ride  comfort  were  evaluated  by  having  two  operators  ride  in  the 
survey  vehicle. 

The  resulting  performance  is  contained  in  the  table  below. 


Table  2:  Detailed  Performance  Criteria 


Performance  Criteria 

Expected  Performance 
Metric  (pre  demo) 

Performance 
Confirmation  Method 

Actual  (post  demo) 

PRIMARY  CRITERIA  (Performance  Objectives) 

Small  to  moderate 
directional  effect  of 
vehicle  magnetic 
signature  in  unfiltered 
data 

<25  nT 

Survey  a  5.5  acre 
section  of  the  Legacy 
Area  in  two  different 
configurations;  Measure 
the  strength  of  streaks  in 
uncorrected  data. 

20nT  for  15’  data 

12nT  for  17’  data 

Small  directional 
effect  of  vehicle 
magnetic  signature  in 
filtered  data 

<2nT 

Survey  a  5.5  acre 
section  of  the  Legacy 
Area  in  two  different 
configurations;  Measure 
the  strength  of  streaks  in 
corrected  data. 

<  1  nT 

Small  effect  of 
filtering  on  amplitude 
and  size  of  targets 

<15%  difference 
(filtered  vs  unfiltered)  in 
amplitude  and  size 

Perform  de- median 
filter  in  Oasis;  vary  the 
time  window  length  of 
the  filter;  compare  the 
amplitude  of  nine 
objects 

<  2%  amplitude  difference 
<1%  size  difference 
using  a  nominal  six-second 
window  and  averaging 
across  nine  objects 

Small  pitch  and  roll 
effect  of  vehicle 
magnetic  signature 

<  2  nT  over  smooth 

areas 

<  5  nT  over  rough  areas 

Statistical  comparison 
of  data  over  rough  and 
smooth  areas 

Smooth  <  InT 

Rough  < InT 
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SECONDARY  CRITERIA  (Performance  Objectives) 

Vehicle  Ride 

Acceptable 

Having  operator  drive 
the  vehicle  during  a 
representative  survey 

Acceptable 

Vehicle  Terrain- 
Handling  Capability 

Acceptable 

Having  operator  drive 
the  vehicle  during  a 
representative  survey 

Acceptable 

7.3  Data  Analysis,  Interpretation  and  Evaluation 
7.3.1  Bulk  Vehicle  Signature  -  Octant  Tests 

As  described  above,  octant  data  (fifteen  seconds  of  static  data  acquired  with  the  center 
magnetometer  over  the  same  location  on  the  ground,  and  with  the  vehicle  and  towed  platform 
oriented  magnetic  North,  Northeast,  East,  Southeast,  South,  Southwest,  West,  and  Northwest) 
were  acquired  in  four  configurations: 

•  With  the  magnetometers  at  their  current  15’  distance  from  the  vehicle  and  with  the 
vehicle’s  steel  skidplate  attached 

•  With  the  magnetometers  at  their  current  15’  distance  from  the  vehicle  and  with  the 
vehicle’s  steel  skidplate  removed 

•  With  the  magnetometers  moved  back  an  additional  two  feet  and  with  the  vehicle’s  steel 
skidplate  attached 

•  With  the  magnetometers  moved  back  an  additional  two  feet  and  with  the  vehicle’s  steel 
skidplate  removed 

For  each  of  these  tests,  data  were  collected  both  with  the  engine  off  and  with  the  engine  running. 
When  the  engine  is  running,  a  vehicle’s  alternator  or  generator  usually  appears  as  a  directionally- 
dependent  offset;  by  collecting  data  both  ways  we  can  see  the  magnitude  of  the  effect.  Also,  for 
each  test,  the  GPS  antenna  was  removed  from  its  perch  above  the  center  magnetometer  to  be 
certain  that  it  was  not  effecting  the  measurement  of  the  vehicle’s  signature. 

The  figures  below  show  the  octant  data.  All  of  the  octant  plots  below  are  displayed  to  a  30nT 
vertical  scale.  A  reference  magnetometer  was  placed  just  South  of  the  octant  measurement  area, 
and  was  used  to  measure  the  Earth’s  ambient  field;  this  value  was  then  subtracted  from  the  octant 
values  during  processing  to  remove  drift  and  produce  curves  with  the  same  absolute  reference. 
The  reference  magnetometer  malfunctioned  the  last  afternoon  of  the  demonstration,  so  the  last 
set  of  curves  -  the  17’  data  with  the  skidplate  installed  -  had  no  absolute  reference  data.  For  this 
last  set,  the  data  were  manually  shifted  to  place  the  peak  value  of  the  curve  in  a  similar  location 
to  the  other  curves  (just  above  zero). 
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Orientation,  Degrees 


Figure  5:  Octant  Test,  Mag  15’  Behind  Vehicle,  Skidplate  Removed 


- engine  off 

—  engine  on 


Figure  6:  Octant  Test,  Mag  15’  Behind  vehicle,  Skidplate  Installed 
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Octant  Test,  Mag  17'  Behind  Vehicle,  Skidplate  Removed 


- engine  off 

—  engine  on 


Figure  7:  Octant  Test,  Mag  17’  Behind  Vehicle,  Skidplate  Removed 


Octant  Test,  Mag  17'  Behind  Vehicle,  Skidplate  Installed 


- engine  off 

—  engine  on 


Figure  8:  Octant  Test,  Mag  17’  Behind  Vehicle,  Skidplate  Installed 
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These  data  show  several  things.  Firstly,  the  maximum  and  minimum  values  do  not  occur  in 
opposing  directions.  That  is,  the  maximum  value  occurs  at  the  Southern  orientation,  but  the 
minimum  does  not  occur  at  the  same  orientation;  sometimes  the  minimum  occurs  at  the 
Southwest  orientation,  sometimes  East,  and  sometimes  West.  One  could  argue  that,  in  presenting 
these  data,  one  should  present  the  peak-to-peak  difference  as  the  “true”  measure  of  the  bulk 
signature.  However,  since  real-world  surveying  virtually  always  occurs  in  opposing  sets  of 
orientations  (e.g.,  North  and  South,  never  North  and  West),  representing  bulk  signature  as  the 
difference  between  readings  acquired  in  North  and  South  orientations  remains  a  useful  tool. 
Second,  as  expected,  while  the  vehicle’s  engine  is  running,  its  alternator  is  adding  a 
directionally-dependent  component  to  the  vehicle’s  signature,  contributing  perhaps  5nT  to  the 
bulk  signature.  Third,  the  presence  of  the  welded  steel  skid  plate  is  having  a  definite  effect  on  the 
shape  of  the  octant  data  but  does  not  appear  to  be  grossly  effecting  the  magnitude  of  the  bulk 
signature.  Note  that  a  signature  removal  technique,  whether  it  is  table-driven  (looking  up  a  pre¬ 
calculated  offset  for  a  given  orientation  in  a  table  and  subtracting  it  off)  or  filter-based,  must  be 
able  to  handle  these  situations  and  remove  the  signature  regardless  of  direction. 

We  summarize  these  results  (distance  from  vehicle  to  sensor  array;  skid  plate  on  or  off;  engine 
on  or  off)  in  the  table  below. 


Table  3:  Octant  Test  Summary 


Vehicle-Sensor 

Separation 

(feet) 

Skid  Plate 

Engine 

Running 

Bulk 

Signature 

(South- 

North, 

nT) 

15 

Removed 

No 

12.0 

15 

Removed 

Yes 

15.0 

15 

Installed 

No 

17.7 

15 

Installed 

Yes 

22.1 

17 

Removed 

No 

7.5 

17 

Removed 

Yes 

10.3 

17 

Installed 

No 

11.5 

17 

Installed 

Yes 

13.9 

From  here,  we  can  see  that,  worst  case,  at  15’,  even  with  the  skid  plate  installed  and  the  engine 
running,  the  bulk  signature  is  no  larger  than  25nT.  The  best-case  signature  measured  in  this  way 
was  about  8nT  at  17’  with  the  skid  plate  installed  and  the  engine  off. 

Below  we  compare  these  to  data  acquired  dynamically. 

7.3.2  Bulk  Vehicle  Signature  -  Geophysical  Survey  of  Legacy  Area 

As  per  the  demonstration  plan,  we  surveyed  a  section  of  the  Legacy  Area  in  two  configurations: 

•  An  aggressive  configuration  with  the  magnetometers  at  their  current  15’  distance  from 
the  vehicle  and  with  the  vehicle’s  steel  skidplate  attached 

•  A  conservative  configuration  with  the  magnetometers  moved  back  an  additional  two  feet 
and  with  the  vehicle’s  steel  skidplate  removed 
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We  covered  5.5  acres  in  the  North  section  of  the  APG  site.  This  was  a  somewhat  smaller  area 
than  planned.  This  was  because  the  portion  of  the  Legacy  Area  in  the  Southern  section  of  the 
Standardized  Test  Site  -  the  part  to  the  South  of  the  Indirect  Fire  area  -  was  not  always  available 
due  to  live  weapons  testing  at  a  section  of  Aberdeen  Proving  Grounds  just  south  of  the 
Standardized  Test  Site.  However  it  was  more  than  enough  area  to  collect  sufficient  data  for  the 
Performance  Objectives. 

Below  we  show: 

•  Data  acquired  in  2006  with  the  Chenowth  vehicle 

o  with  no  corrections 
o  with  a  bidirectional  correction  applied 
o  with  a  nominal  six-second  median  filter  applied 

•  Data  acquired  in  2009  with  the  COTS  vehicle. 

o  with  no  corrections 
o  with  a  bidirectional  correction  applied 
o  with  a  nominal  six-second  median  filter  applied 

o  with  the  median  filter  window  length  varied  to  measure  the  effect  on  target 
amplitude  and  spatial  extent 

2006  Data  with  Chenowth  Vehicle 

The  scale  for  all  image  displays  below  is  a  relatively  tight  +-  25nT. 

First  we  examine  uncorrected  data  acquired  with  VSEMS  in  2006  using  the  Chenowth  aluminum 
buggy. 
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402700  402750  402800  402850  402900 


APG  Standardized  UXO  Test  Site 
Surveyed  with  Chenowth  Buggy,  2006 

_ 15"  Mag-To- Vehicle  Separation _ 

Figure  9:  APG  Site  Surveyed  in  2006  with  Chenowth  Buggy,  Mags  15’  Back,  Uncorrected 

On  the  western  edge  of  this  data  set,  a  particular  set  of  traverses  clearly  has  a  data-specific 
streaking  issue.  Field  notes  from  2006  indicate  that  a  GPS  antenna  cable  directly  over  the  center 
magnetometer  had  failed,  and  was  replaced  with  a  cable  that  apparently  had  a  ferrous  component 
on  the  coupler,  causing  the  visually  apparent  streaking  in  one  set  of  data  before  it  was 
discovered.  For  simplicity,  we  first  exclude  these  data  and  examine  a  subset  of  the  data  set 
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overlaying  the  northern  section  of  the  Legacy  area.  The  resulting  uncorrected  data  below  are  of 
high  quality  and  show  only  light  streaking  at  this  relatively  tight  +-25nT  display  scale,  as  one 
would  expect  of  data  acquired  using  an  exotic  aluminum  vehicle. 


402750  402800 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  Chenowth  Buggy,  2006 
15"  Mag-To-Vehicle  Separation 

Figure  10:  APG  Site  Legacy  Area  Section  Surveyed  in  2006  with  Chenowth  Buggy,  Mags  15’  Back, 

Uncorrected 

Next  we  apply  a  bidirectional  correction  to  these  data.  This  technique,  originally  developed  by 
GEO-CENTERS  to  remove  light  streaks  from  STOLS  data,  is  not  a  filtering  technique;  it  is  an 
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offset  subtraction  technique.  It  involves  selecting  an  area  of  interest  where  there  are  no 
anomalies,  separating  the  survey  data  lines  into  sets  of  opposite  directions  (e.g.,  for  the  Aberdeen 
data,  North-going  traverses  and  South-going  traverses),  averaging  the  portions  of  data  passing 
through  the  area  of  interest  to  calculate  a  separate  offset  for  each  magnetometer  and  for  each 
direction,  and  then  subtracting  this  offset  from  each  magnetometer  for  each  direction.  The 
individual  magnetometer  offsets  are  shown  in  the  table  below.  For  these  and  all  other 
bidirectional  corrections  referred  to  in  this  document,  the  “_ns”  data  sets  were  acquired  with 
traverses  going  from  North  to  South;  the  “_sn”  data  sets  were  acquired  with  traverses  going  from 
South  to  North. 


Table  4:  Bidirectional  Corrections  for  Unfiltered  2006  Data  Acquired  with  Chenowth  Buggy  at  15’ 


magi 

mag2 

mag3 

mag4 

mag5 

average 

bulk 

signature 

2006  Unfiltered 

apg06_sn 

1.8 

0.5 

0.3 

-2.3 

-0.6 

-0.1 

apg06_ns 

-3.8 

-1.9 

-3.2 

-6.1 

-4.4 

-3.9 

-3.8 

From  the  table,  we  see  that  the  average  offset  of  the  South  to  North  traverses  was  about  0  nT, 
and  the  average  bias  of  the  North  to  South  traverses  was  about  -4  nT.  Thus  the  measured  bulk 
signature  of  the  Chenowth  vehicle  in  these  data  has  a  magnitude  of  about  4  nT.  This  signature 
produced  the  minor  streaking  in  the  image  above.  The  streaking  has  been  completely  removed  in 
the  image  below  using  the  bidirectional  correction.  This  shows  that  data  acquired  using  an  exotic 
vehicle  still  has  some  degree  of  vehicle  signature  in  it,  and  that  the  data  quality  are  improved  by 
an  offset-based  technique. 
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(meters) 

Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  Chenowth  Buggy,  2006 
_ 15”  Mag-To- Vehicle  Separation _ 

Corrected  Bidirectionally 


Figure  11:  APG  Site  Legacy  Area  Section  Surveyed  in  2006  with  Chenowth  Buggy,  Mags  15’  Back,  with 

Bidirectional  Correction 


Next  we  include  the  streaky  section  to  the  west  of  the  above  data  set,  and  apply  the  same 
bidirectional  correction  to  it.  Not  surprisingly,  it  doesn’t  work,  because  the  level  of  streaking  is 
different  from  the  rest  of  the  data  set  -  one  would  need  to  attempt  to  correct  it  individually.  This 


29 


shows  a  major  limitation  of  the  bidirectional  technique  -  that  it  requires  the  individual  offsets  to 
be  unchanging  for  each  direction. 


402800  - 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  Chenowth  Buggy,  2006 
_ 15"  Mag-To- Vehicle  Separation _ 


402750 


402800 


Figure  12:  APG  Site  Legacy  Area  Plus  Adjacent  Section  Surveyed  in  2006  with  Chenowth  Buggy,  Mags  15’ 

Back,  with  Bidirectional  Correction 


Lastly,  we  correct  the  entire  data  set  with  a  de-median  filter  with  a  nominal  six-second  window. 
This  removes  all  streaking,  including  the  streaks  in  the  section  of  data  with  the  errant  cable,  and 
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in  the  Southern  section  which  is  acquired  at  a  completely  distinct  set  of  directions.  This  shows 
the  utility  of  the  filtering  method  over  the  offset-based  method. 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  Chenowth  Buggy,  2006 
15'  Mag-To- Vehicle  Separation 


6  Second  De-median  Filter 


Figure  13:  APG  Site  Legacy  Area  Plus  Adjacent  Section  Surveyed  in  2006  with  Chenowth  Buggy,  Mags  15’ 

Back,  with  De-median  Correction 


Visually,  there  are  no  streaks  in  these  data,  but  we  can  apply  the  same  method  of  measurement  as 
we  did  above:  dividing  the  data  into  South-to-North  traverses  and  North-to-South  traverses, 
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evaluating  the  readings  per  sensor  in  a  quiet  area  of  interest,  and  averaging  the  results.  We  do  not 
apply  the  bidirectional  corrections;  we  merely  generate  them  to  populate  the  table  below,  which 
shows  that  the  bulk  signature  remaining  after  filtering  is  essentially  zero. 


Table  5:  Bidirectional  Corrections  for  Filtered  2006  Data  Acquired  with  Chenowth  Buggy  at  15’ 


magi 

mag2 

mag3 

mag4 

mag5 

average 

bulk 

signature 

2006  Filtered 

apg06_sn 

0.1 

0.1 

0.1 

0.0 

0.0 

0.1 

apg06_ns 

0.0 

0.1 

0.1 

0.1 

0.1 

0.1 

0.0 

2009  Data  with  COTS  Vehicle  -  Aggressive  Configuration  (magnetometers  15’  back, 
skidplate  installed) 

Next  we  look  at  data  acquired  in  2009  over  the  Legacy  Area  with  the  COTS  vehicle  in  its  more 
aggressive  configuration  (with  the  magnetometers  15’  behind  the  vehicle  and  with  the  steel  skid 
plate  installed).  The  streaks  from  the  directionally-varying  signature  in  adjacent  traverses 
acquired  180  degrees  apart  are  plainly  visible.  Unlike  the  uncorrected  2006  Chenowth  data 
which  had  only  minor  streaking,  the  quality  of  the  uncorrected  data  acquired  with  the  COTS 
vehicle  is  clearly  impaired  without  some  sort  of  streak  removal. 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  wth  COTS  Vehicle,  2009 

_ 15"  Mag-To- Vehicle  Separation,  Plate  Attached _ 

Figure  14:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15'  Back,  Skidplate  Attached,  Uncorrected 

Next  we  apply  the  same  kind  of  bidirectional  streak  correction  described  above  to  the  2009 
COTS  vehicle  data.  Unlike  the  2006  data  which  were  acquired  in  discrete  lines  (pausing  data 
acquisition  while  the  vehicle  turned  around),  the  2009  data  was  acquired  in  one  long  continuous 
line,  with  the  turnarounds  intentionally  kept  in  the  data  set  so  we  could  see  the  effect  of  the 
vehicle.  To  perform  the  bidirectional  correction,  we  first  need  to  run  a  program  that  chops  up  the 
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single  long  data  line  into  distinct  lines  in  each  of  the  primary  directions  (Oasis  will  do  this  as 
well,  through  the  Database  Tools  I  Line  Tools  I  Split  on  XY  option).  Then  the  correction  is 
performed  as  described  above.  The  offset  values  are  shown  in  the  table  below. 


Table  6:  Bidirectional  Corrections  for  Unfiltered  2009  COTS  Vehicle  Data  Acquired  at  15’  with  Skidplate 


magi 

mag2 

mag3 

mag4 

mag5 

average 

bulk 

signature 

2009  Aggressive 

Unfiltered 

apg09_1 5p_sn 

-22.8 

-22.6 

-24.3 

-32.5 

-28.0 

-26.1 

apg09_1 5p_ns 

-7.3 

-0.8 

-3.6 

-11.2 

-7.3 

-6.0 

20.0 

From  here  we  see  that,  using  the  COTS  vehicle  in  its  more  aggressive  configuration,  the  average 
bias  of  the  South-to-North  traverses  is  about  -26  nT;  the  average  bias  of  the  North-to-South 
traverses  is  about  -6  nT.  This  yields  a  bulk  signature  of  about  20  nT,  as  compared  with  4  nT  with 
the  Chenowth  vehicle.  With  the  individual  biases  subtracted  off,  the  resulting  data  in  the  image 
below  are  streak-free.  This  provides  confirmation  that  high  quality  magnetometer  data  can  be 
acquired  using  a  COTS  vehicle  and  applying  a  simple  bidirectional  offset-based  correction. 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
_ 15"  Mag-To- Vehicle  Separation,  Plate  Attached _ 

_ Corrected  Bidirectionally _ 

Figure  15:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15'  Back,  Skidplate  Attached,  with  Bidirectional 

Correction 

The  usefulness  of  the  bidirectional  correction,  however,  has  limits.  Firstly,  it  requires  data  to  be 
collected  in  sets  of  uniform  direction.  This  is  not  a  problem  for  grid  surveys  in  obstructed  areas, 
but  in  real-world  production  geophysics,  sites  are  sometimes  oddly-shaped,  have  obstacles  that 
must  be  driven  around,  or  require  patchwork  fill-ins  that  make  the  bidirectional  technique 
cumbersome  or  useless.  In  addition,  experience  has  shown  that,  even  on  large  regularly-shaped 
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sites  where  physical  tire  tracks  on  the  ground  are  the  main  source  of  guidance,  the  survey 
direction  will  not  remain  invariant  and  will  instead  begin  to  drift  as  tire  tracks  are  imperfectly 
followed.  These  new  tracks  will  be  at  a  slightly  different  angle,  and  one  set  of  biases  will  not  be 
able  to  take  out  all  of  the  streaks.  This  raises  the  second  issue.  Because  the  bidirectional 
correction  is  an  offset  subtraction  technique,  not  a  filtering  technique,  it  assumes  that  the  offset  is 
not  changing  at  all,  and  thus  fails  if  the  offset  changes,  either  due  to  physical  changes  on  the 
platform  (such  as  the  connector  issue  in  2006)  or  if  the  survey  direction  changes.  A  substantial 
advantage,  then,  to  filtering  techniques  is  that  they  can  not  only  remove  a  directional  offset,  they 
can  remove  a  changing  directional  offset.  Of  course,  they  can  also  remove  long-wavelength 
effects  from  geology  and  man-made  structures.  For  example,  although  the  bidirectional 
correction  has  nicely  removed  the  streaks  from  the  above  image,  it  has  done  nothing  to  remove 
the  positive  bloom  in  the  northeast  corner  of  the  image  that  is  likely  created  by  a  combination  of 
site  geology,  proximity  to  a  railroad  track,  and  an  overhead  power  line.  Below,  we  apply  a  six- 
second  de-median  filter  to  the  data.  The  degree  to  which  it  removes  not  only  the  vehicle 
signature  but  the  bloom  in  the  northeast  corner  is  visually  apparent. 

In  addition  to  the  image  being  visually  clear  of  streaks,  we  perform  the  bidirectional  analysis  on 
the  filtered  data  to  generate  the  corrections  in  the  table  below,  which  shows  that,  after  filtering, 
the  remaining  bulk  signature  is  essentially  zero. 


Table  7:  Bidirectional  Corrections  for  Filtered  2009  COTS  Vehicle  Data  Acquired  at  15’  with  Skidplate 


magi 

mag2 

mag3 

mag4 

mag5 

average 

bulk 

signature 

2009  Aggressive 

Filtered 

apg09_1 5p_sn 

0.1 

0.0 

0.0 

0.1 

0.1 

0.0 

apg09_1 5p_ns 

0.1 

0.1 

0.1 

0.1 

0.0 

0.1 

0.1 
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402660  402680  402700  402720  402740  402760  402780  402800 


(meters) 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

Six  Second  De-median  Filter 


Figure  16:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  De-median 

Filter  Correction 


2009  Data  with  COTS  Vehicle  -  Conservative  Configuration  (magnetometers  17’  back, 
skidplate  removed) 

Next  we  examine  the  2009  COTS  vehicle  data  in  its  more  conservative  configuration  -  with  the 
magnetometers  moved  back  to  17’,  and  with  the  steel  skidplate  removed.  The  image  below 
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shows  some  light  streaking,  but  the  streaks  are  of  a  lower  magnitude  than  those  in  the  more 
aggressive  configuration  above. 


402675  402700  402725  402750  402775  402800 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
17'  Mag-To- Vehicle  Separation,  Plate  Removed 


Figure  17:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  17'  Back,  Skidplate  Removed,  Uncorrected 

We  divided  the  data  into  directionally  uniform  sets,  applied  a  bidirectional  correction,  and 
generated  the  values  in  the  table  below  showing  a  bulk  signature  of  approximately  12  nT  in  this 
configuration. 
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Table  8:  Bidirectional  Corrections  for  Unfiltered  2009  COTS  Vehicle  Data  Acquired  at  17’  without  Skidplate 


magi 

mag2 

mag3 

mag4 

mag5 

average 

bulk 

signature 

2009  Conservative  Unfiltered 

apg09_17n_sn 

-10.2 

-12.1 

-9.2 

-17.4 

-12.3 

-12.3 

apg09_17n_ns 

0.9 

1.4 

2.7 

-6.0 

-1.1 

-0.4 

11.9 

Data  with  the  bidirectional  correction  applied  are  shown  below. 
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402675  402700  402725  402750  402775  402800 


(meters) 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
17  Mag -To- Vehicle  Separation,  Plate  Removed 

Bidirectional  Correction 


Figure  18:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  17’  Back,  Skidplate  Removed,  with  Bidirectional 

Correction 
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The  bidirectionally-corrected  data  shown  above  exhibit  some  light  streaking  visible  near  the 
center  of  the  image.  We  attribute  these  streaks  to  the  reliance  of  the  correction  on  rigorously 
uniform  direction  in  both  the  survey  data  and  in  the  data  used  to  generate  the  correction.  Below 
we  apply  the  standard  six-second  de-median  filter  instead;  it  removes  the  light  streaking  in  the 
middle,  and  also  removes  most  of  the  cultural  and  geological  bloom  in  the  northeast  corner. 
Below  we  show  the  bulk  signature  data,  again  showing  that  after  filtering,  the  bulk  signature  is 
essentially  zero. 


Table  9:  Bidirectional  Corrections  for  Filtered  2009  COTS  Vehicle  Data  Acquired  at  17’  without  Skidplate 


magi 

mag2 

mag3 

mag4 

mag5 

average 

bulk 

signature 

2009  Conservative 

Filtered 

apg09_17n_sn 

0.0 

0.0 

0.0 

0.0 

0.1 

0.0 

apg09_17n_ns 

0.0 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
17  Mag -To- Vehicle  Separation,  Plate  Removed 

6  Second  De-median  Filter 


Figure  19:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  17’  Back,  Skidplate  Removed,  with  De-median 

Filter  Correction 


7.3.3  Bulk  Signature  Summary 

In  the  table  below,  we  summarize  the  bulk  signature  of  the  octant  data  and  the  survey  data 
(uncorrected  and  de-median  filtered  with  a  nominal  six-second  window)  for  the  configurations 
we  have  described  above.  We  position  the  octant  and  survey  data  measurements  from  like 
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configurations  in  adjacent  lines  so  their  consistency  can  be  seen.  As  per  the  above  tables,  the 
remaining  bulk  signature  after  filtering  is  essentially  zero. 


Table  10:  Bulk  Signature  Summary 


Vehicle-Sensor  Separation  (feet) 

Octant  or 
Survey 

Bulk  Signature, 
Unfiltered  (nT) 

Bulk  Signature, 
Filtered  (nT) 

2006,  15',  Chenowth 

Survey 

3.8 

0.0 

2009,  15',  COTS,  Skidplate  Removed 

Octant 

15.0 

2009,  15',  COTS,  Skidplate  Installed 

Octant 

22.1 

2009,  15',  COTS,  Skidplate  Installed 

Survey 

20.0 

0.1 

2009,  17',  COTS,  Skidplate  Removed 

Octant 

10.3 

2009,  17',  COTS,  Skidplate  Removed 

Survey 

11.9 

0.1 

2009,  17',  COTS,  Skidplate  Installed 

Octant 

7.9 

7.3.4  Effect  of  De-median  Filter  Window  Length 

We  know  that  the  streaks  in  the  data  are  caused  by  the  vehicle’s  magnetic  signature  changing 
with  the  vehicle’s  direction  relative  to  North.  But  irrespective  of  the  cause,  as  long  as  the 
magnitude  of  the  streaks  are  changing  slowly  relative  to  the  objects  of  interest  in  the  data,  the 
streaks  can  easily  be  filtered  out  using  a  de-median  filter.  A  de-median  filter  is  a  high-pass  filter 
that  calculates  the  median  (middle)  value  for  a  window  of  data  points,  then  subtract  this  median 
value  from  the  value  at  the  center  of  the  window.  This  acts  like  a  dynamic  background 
subtraction,  removing  not  only  changing  remnant  vehicle  signature  but  also  diurnal  drift,  long- 
wavelength  geology,  and  small  intermagnetometer  offsets.  We  did  not  pioneer  the  use  of  de¬ 
median  filters  for  towed  magnetometer  array  data  (we  believe  they  were  first  applied  to  vehicular 
MTADS  data  by  Dan  Steinhurst  and  Dave  Wright  in  the  1990s),  but  we  have  employed  a  de¬ 
median  filter  since  VSEMS  was  first  deployed  in  2002.  Above  we  demonstrated  the  successful 
application  of  a  de-median  filter  with  a  six-second  window  to  the  APG  2006  data  acquired  with 
the  Chenowth  vehicle,  and  to  the  2009  data  with  the  COTS  vehicle.  A  de-median  filter  has  the 
additional  utility  that  it  is  implemented  in  Geosoft  Oasis.  It  is  accessible  in  two  places:  through 
the  UxProcess  module  (under  “Instrument  Drift  Correction,”  where  the  window  size  is  specified 
in  seconds),  and  through  the  Database  Tools  I  Filters  menu  (under  “Rolling  Statistics,”  where  the 
window  size  is  specified  in  fiducials;  note  that  the  program  calculates  the  median,  and  one  must 
manually  set  up  a  separate  channel  and  subtract  the  median  from  the  original  data  to  calculate  the 
de-median).  These  options  allow  any  Oasis  user  to  easily  de-median  filter  magnetometer  data 
acquired  with  a  COTS  vehicle. 

As  with  any  filtering  technique,  successful  use  of  a  de-median  filter  requires  that  the  signals  of 
interest  (in  our  case,  buried  ferrous  objects  of  a  size  and  shape  similar  to  MEC)  have  a 
wavelength  that  is  distinct  from  that  which  one  wishes  to  filter  out.  The  length  of  the  window 
used  for  a  de-median  filter  is  a  tradeoff  between  preservation  of  signal  and  responsiveness  to  the 
changing  directional  effects  of  the  vehicle  signature.  A  long  window  will  preserve  the  size  and 
amplitude  of  MEC  signals  but  will  not  be  responsive  when  the  vehicle  changes  direction  quickly 
(such  as  when  the  vehicle  turns).  Conversely,  using  a  very  short  window  will  produce  a  filter  that 
is  highly  responsive  to  fast  changes  in  vehicle  direction  but  it  may  adversely  affect  MEC  signal. 
Historically  we  have  defaulted  to  a  median  filter  length  of  six  seconds,  a  value  which  was  chosen 
by  Edisonian  principles  (trial  and  error).  For  this  analysis,  we  track  the  magnitude  and  spatial 
extent  of  nine  anomalies  as  we  alter  the  filter’s  window  length  from  1  to  15  seconds. 
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For  completeness,  we  include  all  15  median- filtered  images  below,  starting  with  a  one-second 
window,  and  ending  with  a  15  second  window.  Like  the  prior  images,  these  images  are  at  a 
display  scale  of  +-25  nT.  It  is  visually  apparent  that  the  one  and  two  second  windows  are 
affecting  the  shape  of  anomalies  (below,  we  quantify  this  by  measuring  the  anomaly  full  width  at 
half  max).  It  is  also  clear  that  a  certain  amount  of  streaking  has  crept  back  into  the  images  at  the 
longer  windows  due  to  the  unresponsiveness  of  the  filter  to  changing  vehicle  direction. 
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402660  402680 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To-Vehicle  Separation,  Plate  Installed 

1  Second  De-median  Filter 


Figure  20:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  1  Second  De- 

median  Filter  Correction 
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402660  402680  402700  402720  402740  402760  402780  402800 


402660  402680  402700  402720  402740  402760  402780  Sca/e^§32  7586 


I  I  I 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 


2  Second  De-median  Filter 


Figure  21:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  2  Second  De- 

median  Filter  Correction 
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402660  402680  402700  402720  402740  402760  402780  402800 


(meters) 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To-Vehicle  Separation,  Plate  Installed 

3  Second  De-median  Filter 


Figure  22:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  3  Second  De- 

median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

4  Second  De-median  Filter 


Figure  23:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  4  Second  De- 

median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

5  Second  De-median  Filter 


Figure  24:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  5  Second  De- 

median  Filter  Correction 
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402660  402680  402700  402720  402740  402760  402780  402800 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

6  Second  De-median  Filter 


Figure  25:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  6  Second  De- 

median  Filter  Correction 
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402660  402680  402700  402720  402740  402760  402780  402800 


(meters) 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

7  Second  De-median  Filter 


Figure  26:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  7  Second  De- 

median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 


Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

8  Second  De-median  Filter 


Figure  27:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  8  Second  De- 

median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

9  Second  De-median  Filter 


Figure  28:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  9  Second  De- 

median  Filter  Correction 
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(meters) 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 


10  Second  De-median  Filter 


Figure  29:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  10  Second 

De-median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

11  Second  De-median  Filter 


Figure  30:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  11  Second 

De-median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
1 5' Mag -To -Vehicle Separation,  Plate  Installed 

12  Second  De-median  Filter 


Figure  31:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  12  Second 

De-median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

13  Second  De-median  Filter 


Figure  32:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  13  Second 

De-median  Filter  Correction 
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Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 


14  Second  De-median  Filter 


Figure  33:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  14  Second 

De-median  Filter  Correction 
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Section  of  APG  Standardized  UXOTestSite 
Surveyed  with  COTS  Vehicle,  2009 
15'Mag-To-VehicleSeparation,  Plate  Installed 

15  Second  De-median  Filter 


Figure  34:  Legacy  Area  Surveyed  with  COTS  Vehicle,  Mags  15’  Back,  Skidplate  Installed,  with  15  Second 

De-median  Filter  Correction 


To  quantify  the  effect  of  the  filtering,  we  selected  nine  targets  of  a  variety  of  sizes  (spatial 
anomaly  extent)  from  the  data  set.  For  each  target,  for  each  of  the  15  median  filter  window  filter 
lengths,  we  extracted  the  best  time-series  profile  over  that  target,  and  tabulated  the  target’s 
amplitude  and  the  full  width  at  half  the  maximum  value.  These  are  shown  in  both  tabular  and 
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graphical  form  below.  When  looking  at  the  peak  amplitude,  the  graphs  level  out  for  most  of  the 
nine  objects  at  approximately  three  seconds.  The  exception  is  object  #2,  which  is  the  object 
whose  anomaly  has  the  largest  spatial  extent;  the  curve  for  this  object  flattens  after 
approximately  six  seconds.  To  evaluate  compliance  with  the  performance  objective  of  not 
altering  the  anomaly  by  more  than  15%,  we  take  the  most  slowly  changing  item  (item2);  its 
strength  is  within  15%  of  its  asymptotic  value  after  four  seconds,  and  within  1%  after  6  seconds. 


Table  11:  Nine  Objects,  Anomaly  Peak  Amplitude  as  a  Function  of  Median  Filter  Window  Length 


iteml 

item2 

item3 

item4 

item5 

item  6 

item  7 

item8 

item9 

1 

458.5 

385.5 

1247.2 

206.1 

3861.4 

369.4 

141.1 

510.2 

322.5 

2 

1009.1 

874.0 

1959.7 

546.4 

6841.7 

914.7 

375.5 

854.7 

801.5 

3 

1386.3 

1376.3 

2146.2 

740.6 

7490.8 

1180.7 

505.7 

969.0 

1091.6 

4 

1579.8 

1682.7 

2197.2 

833.7 

7616.3 

1285.4 

553.6 

1013.8 

1219.4 

5 

1659.9 

1871.5 

2214.7 

881.9 

7647.0 

1328.9 

565.4 

1032.0 

1255.5 

6 

1691.6 

1972.0 

2217.2 

906.0 

7649.7 

1353.2 

568.8 

1041.0 

1272.8 

7 

1707.7 

2022.9 

2218.8 

920.7 

7649.7 

1369.2 

568.8 

1042.6 

1276.7 

8 

1714.6 

2052.0 

2220.4 

925.3 

7649.7 

1378.0 

568.8 

1042.6 

1277.7 

9 

1718.2 

2066.8 

2222.9 

928.0 

7649.7 

1379.6 

568.8 

1042.6 
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Figure  35:  Nine  Objects,  Anomaly  Peak  Amplitude  as  a  Function  of  Median  Filter  Window  Length 
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For  anomaly  size  (full  width  at  half  max),  most  of  the  graphs  visually  level  out  at  approximately 
six  seconds,  with  the  graphs  of  objects  1,  2,  and  3  taking  slightly  longer  to  flatten.  To  evaluate 
the  performance  objective  of  not  altering  the  anomaly  by  more  than  15%,  we  take  the  most 
slowly  changing  item  (item2);  its  size  is  within  15%  of  its  asymptotic  value  after  four  seconds, 
and  within  1%  after  6  seconds. 


Table  12:  Nine  Objects,  Anomaly  Full  Width  at  Half  Max  as  a  Function  of  Median  Filter  Window  Length 


seconds 

iteml 

item2 

item3 

item4 

item5 

item  6 

item  7 

item8 

item9 

1 

0.34 

0.35 

0.39 

0.66 

0.5 

0.46 

0.46 

0.38 

0.52 

2 

0.84 

0.74 

0.53 

1.23 

0.79 

0.86 

0.82 

0.6 

0.84 

3 

1.11 

1.14 

0.56 

1.4 

0.83 

1.02 

0.98 

0.65 

1.03 

4 

1.22 

1.29 

0.6 

1.52 

0.86 

1.1 

1.03 

0.67 

1.11 

5 

1.28 

1.4 

0.6 

1.78 

0.86 

1.12 

1.05 

0.69 

1.13 
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1.3 
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0.6 
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0.86 

1.12 

1.05 

0.71 

1.26 

7 

1.3 
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0.6 
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0.86 

1.14 

1.05 

0.71 

1.26 
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0.6 

1.85 
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0.71 
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Figure  36:  Nine  Objects,  Anomaly  Full  Width  at  Half  Max  as  a  Function  of  Median  Filter  Window  Length 
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To  summarize  the  above  results,  we  take  the  two  tables,  calculate  each  value  as  a  percentage  of 
the  asymptotic  value,  and  average  across  all  nine  objects.  The  six-second  value  nominally  used  is 
highlighted. 

Table  13:  Summary  of  Filter  Effects  on  Amplitude  and  Width 


window 

amplitude 

width 

1 

66.8% 

57.0% 

2 

31.9% 

26.2% 

3 

14.2% 

13.1% 

4 

6.4% 

6.9% 

5 

3.1% 

3.2% 

6 

1.6% 

1.0% 

7 

0.8% 

0.5% 

8 

0.5% 

0.1% 

9 

0.3% 

0.0% 

10 

0.2% 

0.0% 

11 

0.1% 

0.0% 

12 

0.1% 

0.0% 

13 

0.1% 

0.0% 

14 

0.0% 

0.0% 

15 

0.0% 

0.0% 

Together,  these  tables  and  graphs  show  that  both  the  amplitude  and  spatial  extent  of  the  de¬ 
median  filtered  anomalies  meet  the  demonstration’s  metric  (they  are  within  15%  of  the 
asymptotic  value)  when  a  de-median  filter  of  at  least  a  four-second  time  window  is  used.  Further, 
both  the  amplitude  and  size  are  within  2%  of  their  asymptotic  values  when  a  six-second  window 
is  used.  Certainly  using  a  median  filter  in  the  six-second  range  does  not  have  a  deleterious  effect 
on  anomaly  amplitude  and  size. 

As  the  filter  length  increases,  the  question  is  not  whether  there  is  too  much  filtering,  but  whether 
there  is  too  little.  This  is  a  more  subjective  question  than  the  minimum  filter  length,  as  it  involves 
parameters  that  vary  across  the  data  set,  such  as  anomaly  density  and  how  fast  the  vehicle 
changes  direction.  Coincidentally,  the  Chenowth  vehicles  have  a  poor  (large)  turning  radius,  and 
the  COTS  Club  Car  1550XRT  vehicle  has  an  extremely  small  turning  radius.  Thus,  by  using  this 
new  COTS  vehicle,  we  have  introduced  a  vehicle  with  both  a)  a  higher  magnetic  signature  and  b) 
the  ability  to  have  the  amount  of  signature  in  the  data  change  very  quickly.  Further,  in  addition  to 
allowing  the  signature  to  change  very  quickly,  a  small  turning  radius  radically  changes  the  towed 
array  geometry  by  making  it  so  that  the  platform  is  no  longer  directly  behind  the  vehicle;  the 
magnetometers  on  the  side  in  the  direction  of  the  turn  are  now  considerably  closer  to  the  vehicle 
than  those  on  the  other  side.  In  the  two  images  below  we  zoom  in  on  an  area  where  the  vehicle  is 
turning  around  in  as  small  a  space  as  possible.  The  four  black  squares  show  artifacts  of  the  rapid 
turn  that  are  plainly  visible  on  the  inboard  magnetometers.  The  left  image  was  created  using  a 
two-second  de-median  filter,  which  according  to  the  above  analysis  is  too  short  (i.e.,  we  know 
that  it  has  deleterious  effects  on  the  data),  and  the  right  image  uses  a  four-second  filter,  which 
according  to  the  above  analysis  is  barely  adequate.  The  “too  short”  filter  is  removing  some  but 
not  all  of  the  effect,  and  the  “barely  adequate”  filter  is  not  able  to  remove  the  effect. 
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Figure  37:  Section  of  rapid  vehicle  turnaround,  showing  artifacts  on  inboard  mags,  with  2-second  filter  (left) 

and  4-second  filter  (right) 


On  a  standard  vehicular  survey,  when  we  acquire  data,  we  normally  pause  the  system  while  the 
vehicle  is  turning  around  and  do  not  store  the  turnaround  data.  Other  contractors  acquire  the  data 
during  turnarounds  but  chop  it  out  during  processing.  We  acquired  turnaround  data  during  the 
APG  survey  precisely  to  be  able  to  see  this  effect.  It  would  be  unfair  to  tar  the  COTS  vehicle  for 
having  maneuverability  that  the  Chenowth  vehicle  does  not  have.  It  is  important  to  note  that  the 
lesson  here  is  not  that  the  COTS  vehicle  has  too  high  a  signature  to  be  useful  or  that  the  filtering 
is  not  adequate;  it  is  that,  as  the  vehicle’s  signature  increases,  care  must  be  taken  to  restrict  data 
acquisition  and  /or  processing  to  data  whose  direction  is  not  changing  fast  enough  to  cause 
artifacts  that  could  be  processed  and  analyzed  as  false  positives. 


7.3.5  Effect  of  Rough  Terrain 

As  we’ve  said,  although  the  effect  of  the  vehicle  changes  primarily  with  the  vehicle’s  orientation 
relative  to  North,  it  also  changes  as  the  vehicle  and  the  sensor  platform  both  pitch  and  roll  over 
rough  terrain.  The  question  is  whether,  for  a  given  vehicle  and  sensor  separation,  pitch  and  roll 
effects  are  large  enough  to  be  problematic.  In  data  acquired  using  a  variety  of  COTS  vehicles  at  a 
variety  of  sensor  separations  during  the  project’s  testing  at  Fort  Devens,  we  found  that,  when  the 
vehicle  had  too  strong  a  signature  or  the  sensors  were  too  close,  we  saw  this  effect  as  a 
“clouding”  of  the  data  -  low-level  spatially  diffuse  anomalies  that  were  plainly  visible  after  the 
slowly-changing  directional  component  had  been  removed  with  a  de-median  filter.  For  the  APG 
demonstration,  we  wished  to  try  to  correlate  such  possible  artifacts  with  rough  patches  of  terrain 
by  acquiring  GPS  coordinates  over  certain  particularly  rough  sections  of  the  site.  This  was  done 
by  having  a  passenger  take  notes  while  the  driver  was  covering  planned  traverses.  For  example,  a 
sample  note  might  read  “north  end  of  line  5  very  rough”.  After  the  survey  was  completed,  the 
driver  went  back  to  the  notated  areas  and  recorded  their  GPS  coordinates  to  see  if  we  could 
correlate  these  rough  sections  with  visible  artifacts  in  the  filtered  data  such  as  the  “cloudy” 
artifacts  we  saw  with  some  vehicles  at  Devens. 


In  the  two  images  below,  the  rough  sections  at  APG  are  shown  as  the  black  lines.  The  images  of 
six-second  de-median  filtered  data,  like  all  other  images  in  this  document,  are  shown  at  +-  25nT, 
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but  the  point  at  which  green  transitions  to  yellow  has  been  pushed  down  to  zero  in  order  to 
visually  highlight  small  positive  changes  in  the  data. 


Section  of  APG  Standardized  UXO  Test  Site 
Surveyed  with  COTS  Vehicle,  2009 
15’  Mag-To- Vehicle  Separation,  Plate  Installed 

6  Second  De-median  Filter 


Figure  38:  Legacy  area  with  rough  areas  shown  as  black  lines,  color  scale  pushed  to  show  washboard  pattern 
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Figure  39:  Blowup  to  show  washboard  pattern 


We  see  no  such  “cloudy”  artifacts  in  the  more  aggressive  15’  data,  nor  in  the  more  conservative 
17’  data.  What  we  do  see,  however,  is  a  down-track  washboard  artifact,  periodic  and  coherent 
across  all  sensors,  of  approximately  3nT  peak  to  peak.  The  second  image  is  zoomed  in  to  better 
show  the  artifact.  The  washboard  pattern  appears  to  be  uncorrelated  with  the  black  lines 
indicating  rough  areas.  The  pattern  is  almost  certainly  due  to  the  ferrous  beads  in  the  tires  on  the 
buggy.  The  physical  distance  between  the  anomaly  peaks  is  approximately  1.8  meters,  which  is 
almost  exactly  the  circumference  of  the  vehicle’s  1.875’  (0.57m)  high  tires.  The  fact  that  the 
artifact  ebbs  and  flows  somewhat  is  probably  due  to  constructive  interference  between  the 
signals  from  the  steel  beads  on  the  vehicle’s  two  rear  tires.  We  had  demagnetized  the  eight  beads 
on  the  four  tires  when  we  purchased  the  vehicle,  but  that  was  nearly  eighteen  months  before  the 
APG  survey.  It  is  likely  that  they  need  to  be  demagnetized  again.  Note  that  both  MTADS’  and 
VSEMS’  practice  has  been  to  employ  tires  with  their  steel  beads  removed  on  the  survey 
platform,  as  these  are  very  close  to  the  magnetometers,  but  the  tires  on  the  survey  vehicle  must 
bear  the  weight  of  the  vehicle,  so  we  have  been  demagnetizing,  not  removing,  the  beads  in  the 
vehicle’s  tires. 

Even  though  we  do  not  see  any  artifacts  that  we  would  associate  with  pitch  and  roll  over  rough 
terrain,  a  statistical  comparison  was  made  between  data  from  an  anomaly-free  area  in  the  rough 
section  in  the  wet  area  at  the  Northern  edge,  and  data  from  an  anomaly-free  area  that  was  not 


65 


flagged  as  rough.  These  sample  areas  were  analyzed  with  the  statistics  tool  in  Oasis.  The  results 
are  in  the  table  below. 

Table  14:  Statistical  Comparison  of  Data  from  Rough  and  Smooth  Terrain 


min 

max 

mean 

stdev 

smooth 

-1.90 

2.48 

0.12 

0.67 

rough 

-2.14 

2.89 

0.04 

0.81 

7.3.6  Vehicle  Ride 

Despite  the  fact  that  the  APG  site  is  extremely  rugged  and  the  vehicle  does  not  have  rear  wheel 
independent  rear  suspension,  the  vehicle’s  ride  was  perfectly  adequate  for  extended  days  of 
production  geophysics. 

7.3.7  Vehicle  Terrain  Handling  Capability 

The  APG  site  is  very  rough,  with  sections  of  plough  and/or  bulldozer  tracks  at  an  angle  relative 
to  the  North-South  traverses  that  we  ran.  In  addition,  we  surveyed  into  the  muddy  wet  section  at 
the  Western  end  of  the  site.  None  of  these  things  proved  difficult  for  the  COTS  vehicle.  Unlike 
the  Chenowth  vehicle,  the  COTS  vehicle  has  four  wheel  drive,  which  was  a  definite  advantage  in 
these  areas.  Another  advantage  was  turning  radius  -  the  COTS  vehicle  has  a  substantially 
smaller  turning  radius  than  the  Chenowth. 

7.3.8  Conclusions 

The  COTS  vehicle  met  all  of  its  metrics  for  a  successful  demonstration.  As  such,  we  have  shown 
that,  for  production  geophysics,  a  carefully- selected  commercial-off-the-shelf  vehicle  can  be 
used  instead  of  an  exotic  custom-fabricated  aluminum  buggy,  and  the  vehicle  signature  can 
either  be  removed  with  the  bidirectional  correction  (essentially  subtracting  a  constant  offset  from 
each  sensor  for  each  of  two  survey  directions),  or  with  a  de- median  filter  of  approximately  six- 
second  length  that  is  already  implemented  in  Geosoft  Oasis  Montaj.  We  saw  nothing  in  the  data 
acquired  with  VSEMS’  standard  15’  magnetometer-to- vehicle  separation  that  forces  us  to  move 
the  magnetometers  further  back  or  remove  any  of  the  COTS  steel  components,  such  as  the 
skidplate,  and  replace  them  with  expensive  custom-fabricated  non-ferrous  pieces. 

Our  result  -  that  the  vehicle’s  20nT  signature  can  be  easily  removed;  that  the  rugged  APG  terrain 
did  not  appear  to  engender  any  artifacts  in  the  data;  and  that,  other  than  placing  reasonable  limits 
on  the  tightness  of  turnarounds,  the  signature  does  not  appear  to  have  deleterious  effects  on  the 
filtered  data  -  is  almost  certainly  applicable  to  other  vehicles  with  a  similar  signature  and  other 
sites  of  similar  roughness. 
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8  Cost  Assessment 


8.1  Cost  Model 

The  cost  model  is  presented  in  the  table  below. 

Table  15:  Cost  Model 


Cost  Element 

Data  Tracked  During  Demonstration 

Estimated  Costs 

Instrument  cost 

Component  costs  and  integration  costs 
•  Engineering  estimates  based  on  current  development 
•  Lifetime  estimate 

Track  consumables  and  repairs 

•  $13,000 

Mobilization  and 
demobilization 

Cost  to  mobilize  to  site 

Derived  from  demonstration  costs 

•  $8,100 

Site  preparation 

No  unique  requirements  encountered 

Instrument  setup 

costs 

Unit:  $  cost  to  set  up  and  calibrate 

Data  Requirements: 

•  Hours  required 

•  Personnel  required 

•  Frequency  required 

•  $275 

•  2  hours 

•  2  people 

•  One-time  setup 

Survey  costs 

Unit:  $  cost  per  hectare 

Data  requirements: 

•  Hours  per  hectare 

•  Personnel  required 

•  $641/hectare 
•  1.8  hours/hectare 

•  2  people 

Detection  data 
processing  costs 

Unit:  $  cost  per  hectare  as  function  of  anomaly  density 

Data  requirements: 

•  Time  required 
•  Personnel  required 

No  detection  performed 

Discrimination  data 
processing 

Unit:  $  per  anomaly 
•  Time  required 
•  Personnel  required 

No  discrimination 
performed 

Instrument  Cost:  This  is  the  cost  of  the  COTS  vehicle.  When  SAIC  purchased  the  Club  Car 
1550XRT  UTV,  it  cost  approximately  $13,000. 

Mob/Demob:  The  cost  to  mobilize  the  vehicle,  the  towed  array,  and  the  two-man  crew  to  APG 
for  the  demonstration  was  approximately  $8,100. 

Site  Preparation:  No  site  preparation  above  that  which  is  necessary  for  any  vehicle-towed  array 
is  required. 

Instrument  Setup  Costs:  Vehicle  towed  arrays  are  typically  shipped  intact.  The  two-man  crew 
has  routinely  set  the  system  up  and  began  collecting  data  within  an  hour  of  arriving  at  a  site. 
Thus  we  estimate  the  cost  as  an  hour  of  two  people. 

Survey  Costs:  The  survey  costs  are  the  same  as  any  towed  array  (or,  more  specifically,  the  same 
as  a  towed  array  using  an  exotic  custom-built  vehicle),  and  these  of  course  depend  on  many 
factors,  including  site  topography,  hours  of  access  to  the  site,  weather,  GPS  problems,  etc.  The 
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daily  rate  for  the  PI  and  an  equipment  operator  in  the  field  was  $2824/day.  We  surveyed  the  5.5 
acre  (2.2  hectares)  sections  of  the  Legacy  Area  in  about  Vi  a  day.  This  was  not  intended  to  be  a 
test  of  efficiency;  survey  lines  were  short,  and  additional  non-survey  tests  such  as  octant  tests 
were  also  run.  Nonetheless,  this  survey  rate  comes  out  to  be  4.4  hectares/day  or  1.8 
hours/hectare,  at  a  cost  of  $641 /hectare. 

Detection  Data  Processing  Costs:  Processing  of  VSEMS  magnetometer  data  is  similar  to 
processing  data  from  any  other  magnetometer  array.  The  data  must  be  de-spiked  and 
background-leveled.  The  use  of  a  COTS  vehicle  requires  filtering  of  the  data  to  remove  the 
directionally- varying  vehicle  signature,  but  as  we  have  shown,  the  same  de-median  filter  that  is 
traditionally  employed  to  remove  diurnal  drift,  long-wavelength  geology,  and  small 
intermagnetometer  offsets  also  removes  vehicle  signature.  These  steps  are  performed  in  Geosoft 
Oasis  Montaj.  VSEMS’  magnetometer  data  requires  the  additional  step  of  notch-filtering  out  the 
instrument-specific  15  Hz  hum  (created  by  the  60  Hz  ambient  electrical  hum  aliasing  at  15  Hz 
because  it  is  sampled  at  75  Hz).  We  usually  perform  this  step  in  our  own  software,  but  Oasis  is 
capable  of  notch- filtering  the  data  as  well.  Because  our  magnetometer  data  is  IPPS-triggered,  it 
never  requires  latency  correction.  The  magnetometer  data  can  be  read  into  Oasis  and  thresholded 
to  generate  a  mag  dig  sheet. 

Discrimination  Data  Processing  Costs:  The  scope  of  this  project,  to  date,  has  not  included 
discrimination  processing. 

8.2  Cost  Drivers 

Any  vehicle  will  require  a  certain  amount  of  time  to  integrate  sensors,  electronics,  GPS,  and  data 
acquisition  components  into,  but  these  are  little  different  on  COTS  versus  custom  vehicles. 

8.3  Cost  Benefit 

The  cost  benefit  comes  from  not  having  to  spend  $100k  on  a  custom-built  vehicle,  and  instead 
being  able  to  spend  closer  to  $  13k  on  a  commercial-off-the-shelf  vehicle. 
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9  Implementation  Issues 

9.1  Environmental  and  Regulatory  Issues 

Because  the  technology  involves  an  already  validated  sensor  for  digital  geophysical  mapping  - 
total  field  magnetometers  -  there  are  no  specific  regulatory  issues  above  those  that  apply  to  all 
DGM  data.  Any  applicable  regulatory  issues  involve  detection  and  discrimination  systems  of  all 
kinds  (i.e.,  how  clean  is  clean,  etc)  and  are  not  specific  to  this  project  or  technology. 

9.2  End-User  Issues 

Because  the  sensors  are  not  only  COTS  but  the  very  sensors  already  well-used  for  MEC  DGM, 
there  should  not  be  serious  impediments  to  use. 

9.3  Relevant  Procurement  Issues 

The  vehicle,  GPS,  and  magnetometers  are  COTS. 

9.4  Availability  of  the  Technology 

Any  contractor  can  purchase  the  specified  COTS  Club  Car  1550XRT  vehicle  (or  a  similar  one) 
and  use,  along  with  COTS  magnetometers  and  a  GPS,  to  tow  a  magnetometer  array  of  their  own 
design. 

9.5  Specialized  Skills  and  Training 

Because  the  sensors  are  not  only  COTS  but  the  very  sensors  already  well-used  for  MEC  DGM, 
training  is  minimal.  However,  due  to  the  number  and  variety  of  sources  of  noise  (of  which 
directionally-varying  magnetic  signature  is  only  one),  background  or  experience  in 
magnetometry  is  necessary  to  build  a  system  that  collects  high-quality  magnetometer  data. 
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10  Points  of  Contact 


POINT  OF 
CONTACT 

Name 

Organization 

Name 

Address 

Phone/F  ax/email 

Role  in  Project 

Rob  Siegel 

SAIC,  Inc 

7  Wells  Ave 

Newton,  MA  02465 

617-618-4662 

Fax:  617-527-7592 
robert.m.sie£el@saic. 

com 

Principle  Investigator, 
MM-0605 

Bob  Selfridge 

US  Army  Corps  of 
Engineers,  Huntsville 
4820  University 

Square 

Huntsville,  AL 
35816-1822 

256-895-1629 
256-895-1737 
bob.  i.  selfridge  @usace 

.armv.mil 

Contracting  Officer’s 
Technical 
Representative 
(COTR),  MM-0605 

Dr.  Herb  Nelson 

ESTCP  Program 

Office 

901  North  Stuart 

Street,  Suite  303 
Arlington,  VA  22203 

(703)  696-2117 

Fax:  (703)  696-2114 
herb.nelson@osd.mil 

ESTCP  MM  Program 
Manager 
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